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Abstract 

o 

A Goursat structure on a manifold of dimension n is a rank two distribu- 
tion V> such that dim P« = i + 2, for i = 0, n — 2, where denotes the 
derived flag of P, which is defined by Z>(°) = P and p( i+1 ) = X>W + [2?W , pW]. 
Goursat structures appeared first in the work of E. von Weber and E. Cartan, 
who have shown that on an open and dense subset they can be converted 
into the so-called Goursat normal form. Later, Goursat structures have been 
studied by Kumpera and Ruiz. Contact structures on three manifolds and 
^ ■ Engel structures on four manifolds are examples of Goursat structures. In the 

paper, we introduce a new invariant for Goursat structures, called the singu- 
larity type, and prove that the growth vector and the abnormal curves of all 
elements of the derived flag are determined by this invariant. Then we show, 
using a generalized version of Backlund's theorem, that abnormal curves of 
all elements of the derived flag do not determine the local equivalence class 
of a Goursat structure if n > 8. We also propose a new proof of a classical 
theorem of Kumpera and Ruiz. All results are illustrated by the n-trailer 
system, which, as we show, turns out to be a universal model for all local 
Goursat structures. 
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Introduction 



Let D be a smooth rank k distribution on a smooth manifold M, that is a map 
that assigns smoothly to each point p in M a linear subspace T>(p) C T p M of 
dimension k. The derived flag of T> is the sequence defined by T>^ = T> and 
= p(<) + f or j > i. A Goursat structure on a manifold M of 

dimension n > 3 is a rank two distribution P such that, for < % < n — 2, the 
elements of its derived flag satisfy dimV^\p) = i+2, for each point p in M. Goursat 
structures were introduced, using the dual language of Pfaffian systems by E. von 
Weber in 1898. The first period of interest in this special class of distributions 
culminated in the work of Cartan and Goursat. A new period was initiated by 
Giaro, Kumpera, and Ruiz at the end of the seventies. A renewal of interest in 
Goursat structures has been observed from the mid of nineties. 

There are at least three reasons explaining those one century long studies. The 
first reason is that any Goursat structure on R n can be locally converted (on an open 
and dense subset, as it was observed only later by Giaro, Kumpera, and Ruiz p0| ) 
into the so-called Goursat normal form, also known as chained form: 

dx n ' ■ Ln dx n - 1 ^ ^ X ^dx 2 ^dx 1 



It seems that Weber [7(| was the first to exhibit this property and, indeed, Gour- 
sat [2~2"| attributes to him this result. In fact, the starting point of Weber's studies 
was the following question: 'When is a given distribution equivalent to Goursat nor- 
mal form?" , which had led him to discover the concept of derived flag. This question 
is very natural because Goursat normal form has a clear geometric interpretation. 
Indeed, let us consider the space J fc (R, R) of fc-jets of maps from R to R. On the one 
hand, a necessary condition for a curve in J fc (R, R) to be a prolongation of a graph 
of a function from R to R is that it is an integral curve of a distribution which, in the 
canonical coordinates of J fc (R, R), is spanned by the Goursat normal form on R fe+2 . 
On the other hand, a necessary and sufficient condition for a diffeomorphism of 
J fc (R, R) to map prolongations of graphs of functions into prolongations of graphs 
of functions is to preserve the distribution spanned by the Goursat normal form 
on R fc+2 . Such diffeomorphisms are called contact transformations [53| of order k 
and have been intensively studied by Backhand ||, and by Lie and Scheffers [£|0 . 

The second reason of interest in Goursat structures has been the classical prob- 
lem, first considered by Monge, of integrating underdetermined differential equations 
(equivalently, Pfaffian systems) without integration. To be more precise, let T> be a 
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rank k distribution on M. The problem (see e.g. pT]] and WW) is to find k smooth 
functions ipi, . . . ,<fk such that any integral curve j(t) of V can be expressed as 
a smooth function of tpi, . . . ,tpk and their time-derivatives along 7(t). The most 
important achievement of the first period of studies on Goursat structures was a 
result of E. Cartan [11|, who showed that a rank two distribution posses the above 
described property if and only if it is transformable into Goursat normal form. 

The third reason of importance of Goursat structures is that they describe the 
nonholonomic constraints of many mechanical systems. For example, the kinemati- 
cal constraints of a passenger car are described by a Goursat structure on M 2 x (S" 1 ) 2 ; 
those of a truck by a Goursat structure on R 2 x (S 1 ) 3 . Moreover, for Goursat struc- 
tures the nonholonomic motion planning problem can be solved explicitly; either by 
transforming them into Goursat normal form, as suggested by Murray and Sastry 
BTJ and |52fl), or by using the concept of flatness, introduced in control the- 



[see e.g. 

ory by Fliess, Levine, Martin, and Rouchon, which is the above described property 
of calculating the trajectories without integration (see e.g. f43| and jl7| ). 



As we said, the second period of studies on Goursat structures began with a work 
of Giaro, Kumpera, and Ruiz p0[ , who observed that there are Goursat structures 
which are not locally equivalent to Goursat normal form. This observation raised the 
problem of classification of Goursat structures and that of finding their invariants, 
and has led Kumpera and Ruiz to write their important paper [32], where they gave 
a complete classification, up to dimension 7, together with a set of general results 
on Goursat structures. 

In the nineties, research on Goursat structures was concentrated around two 

main topics: the classification problem and the nonholonomic motion planning 

problem for mechanical systems described by Goursat structures. Among results 

concerning the classification problem, Murray [|5T]] obtained, using the concept of 

growth vector^], an easily checkable necessary and sufficient condition for a Goursat 

structure to be equivalent to Goursat normal form (his condition simplifies those of 

Libermann |39| and Kumpera and Ruiz p2[). Cheaito and Mormul |12[ corrected 

the classification in dimension 7 (see also |19|) and obtained a complete classification 

in dimension 8. Mormul |59j obtained a complete classification in dimension 9. It 

turns out that this dimension is the highest one in which there is a finite number of 

non-equivalent Goursat structures. Indeed, Cheaito, Mormul, and the authors [13] 

showed that in higher dimensions there are real continuous parameters in the clas- 

1 The Lie flag of a distribution T> is the sequence defined by Vq = V and = X>, + [Do, T>i], 
for i > 1 . The sequence (dim2?i(p))i>o is called the growth vector of T> at p. 
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sification. Note, however, that in each dimension all Goursat structures are finitely 
determined, which implies that there are no functional parameters in the classifica- 
tion. 

Most of the work concerning mechanical control systems described by Goursat 
structures has been motivated by the study of the n-trailer system. It would be 
impossible to give here a complete set of references on this subject. We have thus 
chosen to cite two books |36|, |38|], and to give a few references concerning standard 
control theory problems for the n-trailer system and chained systems. 

The controllability of the n-trailer system has been proved by Laumond both for 
regular |33| and singular [35] configurations. Improved bounds for the nonholonomy 
degree of the n-trailer at singular configurations have been obtained by S0rdalen |5T] , 
Luca and Risler [j4^| , and Jean |j29 |. For regular configurations, an explicit conversion 
of the n-trailer system into chained form has been obtained by S0rdalen |5t| ; for 
singular configurations, an explicit conversion of the n-trailer system into Kumpera- 
Ruiz normal form has been obtained by the authors f55| (see also Section |||). 

Open loop motion planning has been investigated for general nonholonomic sys- 
tems by Brockett 0], Lafferriere and Sussmann | |33|| , and Liu [^TJ. For chained 
systems, these results have been considerably simplified by Murray and Sastry |52 



(using the special properties of chained form) and by Fliess et al. |17| (using the 
concept of flatness). Combined with the conversion of the n-trailer into chained 
form obtained by S0rdalen |K|, they have led to a solution of the nonholonomic 
motion planning problem for the n-trailer system (see e.g. [37], [5j|, and ||67|| ). 

Path tracking of non-abnormal open loop trajectories has been studied by Fliess 
et al. Wh -, Jiang and Nijmeijer [j30]| , and Walsh et al. |59|. Since for chained systems 
constant trajectories (points) are abnormal, the proposed path tracking strategies 
cannot be applied to achieve pointwise stabilization. Indeed, the linearization of a 
chained system around such trajectories is not controllable. The first who observed 
the difficulties of pointwise stabilization for control systems without drift was Brock- 
ett 0. General algorithms for pointwise stabilization of nonholonomic systems can 
be found in the work of Coron fli] , Pomet [j5f|, McCloskey and Murray [Q, and 



Morin and Samson [53]. For chained systems, improved results have been obtained 
by Samson |59|, S0rdalen and Egeland [62], and Teel et al. [36]. These methods have 
been successfully applied to the n-trailer system (see e.g. Samson ]5JJ, S0rdalen and 
Wichlund ||64j| , and the references given there). 
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Our paper is devoted to a study of the geometry of Goursat structures. We will 
introduce a new local invariant for Goursat structures, called the singularity type, 
whose definition is based on the following observation. If P is a Goursat structure 
then each element PW Q f its derived flag contains an involutive subdistribution 
Ci C that has constant corank one in pW and is characteristic for p( l+1 ). The 
singularity type reflects the geometry of incidence between the distributions pW 
and the distributions Cj. Although, as we prove, the singularity type keeps the same 
information about a Goursat structure as the growth vector, that information is 
encoded in the singularity type in a much more systematic and, what is extremely 
important, in a much more geometric way. In particular, the geometric information 
contained in the singularity type enables us to describe completely all abnormal 
curves of all elements of the derived flag. This can be summarized in the following 
Theorem, which is a combination of Theorem [5.6| and Theorem |6.6| , and gives one 
of the main results of the paper. 

Theorem 1 Let P and P be two Goursat structures defined on two manifolds M 
and M , respectively, of dimension n > 3. Fix two points p and p of M and M, 
respectively. The three following conditions are equivalent: 

(i) The singularity type of P at p equals the singularity type of P at p. 

(ii) The growth vector of T> at p equals the growth vector ofT) at p. 

(iii) There exists a diffeomorphism ip, with p = <p(p), between two small enough 
neighborhoods of p and p that transforms the abnormal curves o/P™ into the 
abnormal curves o/pW, for each i > 0. 

An important example that we will use to illustrate our results on Goursat 
structures will be the n-trailer system, that is a mobile robot (unicycle) towing n 
trailers. In the paper we will calculate rigid curves of the n-trailer and give their 
natural mechanical interpretation: they correspond to motions that fix the positions 
of the centers of at least two trailers. We will also show how to transform locally the 
n-trailer system into a Kumpera-Ruiz normal form, and we will prove a surprising 
result stating that any Goursat structure is locally equivalent to the n-trailer system 
around a well chosen point of its configuration space. This result will enable us 



to use for any Goursat structure a deep result of Jean |29| devoted to singular 
configurations of the n-trailer system, in particular we will extend to all Goursat 
structures Jean's formula for the growth vector of the n-trailer system. In our work, 
the singularity type will replace the angles of the n-trailer system that appear in 
Jean's theorem. 
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In the paper, we will propose an inductive procedure of constructing Kumpera- 
Ruiz normal forms of Goursat structures based on two types of prolongations: regu- 
lar and singular. This construction provides a systematic and unifying approach to 
many results of the paper. In particular, it will be used to show that any Goursat 
structure can be brought to a Kumpera-Ruiz normal form; to study generalized con- 
tact transformations, that is transformations which preserve Goursat structures; and 
to define the above mentioned transformations that transform locally the n-trailer 
system into a Kumpera-Ruiz normal form, and, conversely, that convert locally an 
arbitrary Goursat structure into the n-trailer system around a well chosen point of 
its configuration space. 



Recent studies (see [^6[ and |fRj[| ) show that most distributions are determined 
by their abnormal curves. Our complete description of abnormal curves of Goursat 
structures enables us to conclude that this is not the case of Goursat structures. 



Indeed, combining our study with the main theorem of ]13[ leads us to the following 
result: There exist two Goursat structures T> and V> defined at p and p, respectively, 
that are not locally equivalent but for which there exists a diffeomorphism ip, with 
p = (pip), between two small enough neighborhoods of p and p that transforms 
the abnormal curves of T>^> into the abnormal curves of £> (i) , for each % > sec 
Proposition ^6| and Proposition 

The paper is organized as follows. In the first Section we introduce Goursat 
structures, we give some examples in small dimension, and we define Goursat nor- 
mal form. In the second Section, we provide an inductive definition of Kumpera-Ruiz 
normal form. The proposed concept of prolongations enable us to give a new proof 
of the Kumpera-Ruiz theorem, which states that any Goursat structure can be con- 
verted locally into a Kumpera-Ruiz normal form. In the third Section, we introduce 
the n-trailer system and we construct transformations that bring locally the n-trailer 
system into a Kumpera-Ruiz normal form and, conversely, that bring an arbitrary 
Goursat structure into the n-trailer system. In Section four, we introduce our main 
invariant of Goursat structures, namely, the singularity type. We also compute the 
singularity type for Kumpera-Ruiz normal forms and for the n-trailer system. As we 
have said, the singularity type keeps the same information about Goursat structures 
as the growth vector although in both invariants that information is encoded in a 
different way. Section five is devoted to study relations between these two invari- 
ants. In particular, we give a formula to compute the growth vector of an arbitrary 
Goursat structure and another to compute the singularity type using the growth 
vector. In Section six we study abnormal curves of Goursat structures. We give 
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a complete description of absolutely continuous abnormal curves for all elements 
of the derived flag of any Goursat structure. We prove that the whole information 
about all abnormal curves is given by the singularity type. In Section seven we study 
generalized contact transformations, that is transformations which preserve Gour- 
sat structures (also called symmetries) and we give formulas to calculate them from 
first order contract transformations. Those formulas are used to analyze examples 
of Goursat structures that are non-equivalent but that have diffeomorphic abnor- 
mal curves. The paper ends with three Appendices. The first is devoted to a class 
of distributions that, although of rank greater than two, are very close to Goursat 
structures. This class was also studied by E. von Weber. In the second Appendix 
we construct a normal form which we use in our study of rigidity of integral curves 
of Goursat structures. Finally, in the third Appendix, we illustrate through a set 
of figures different configurations of the n-trailer system corresponding to various 
Kumpera-Ruiz normal forms in dimensions 3, 4, 5, and 6. 
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1 Goursat Structures 



1.1 Derived Flag and Goursat Structures 

A rank k distribution T> on a smooth manifold M is a map that assigns smoothly to 
each point pinMa linear subspace T>(p) C T P M of dimension k. In other words, a 
rank distribution is a smooth rank subbundle of the tangent bundle TM. Such 
a field of tangent fc-planes is spanned locally by k pointwise linearly independent 
smooth vector fields fx, ■ ■ ■ , fk on M, which will be denoted by V = (fx, ■ ■ ■ , fk)- 

Two distributions D and T> defined on two manifolds M and M, respectively, 
are equivalent if there exists a smooth diffeomorphism ip between M and M such 
that 

((p*V)ip) =V(p), 

for each point p in M. They are locally equivalent at two points p and p that 
belong to M and M, respectively, if there exists two small enough neighborhoods 
U and U of p and p, respectively, and a diffeomorphism 9? between U and £/ such 
that (pip) = p and (y?*P)(p) = T>(p), for each point p in [/. 

The derived flag of a distribution D is the sequence T>^ C Z^ 1 -* C ■ ■ • defined 
inductively by 

V {0) =V and =P« + [P«,P«], fori>0. (1) 

This sequence should not be confused with the Lie flag (|37|), which will be introduced 
in Section |. 

Definition 1.1 A Goursat structure on a manifold M of dimension n > 3 is a rank 
two distribution X> such that, for < i < n — 2, the elements of its derived flag 
satisfy dimZ?W(p) = j + 2, for each point p in M. 

1.2 Examples of Goursat Structures 

A Goursat structure on a three-manifold is a contact structure, and is locally equiv- 
alent to the distribution spanned by 

( -9_ B 1 g \ 

y 9x3 ' " 9x2 9xi y ' 

which is called Pfaff-Darboux normal form. A Goursat structure on a four-manifold 
is an Engel structure, and is locally equivalent to the distribution spanned by 

( dxl ' X ^~8x~i + X3 9xl + 9x7 ) ' 
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which is called Engel normal form. Observe that, for a generic field of planes D 
on R 3 , we have dimPW(p) = 3, for any point p on an open and dense subset of R 3 ; 
for a generic field of planes T> on IR 4 , we have dim.T>^ 1 \p) = 3 and diraT)^ (p) = 4, 
for any point p on an open and dense subset of R 4 . Therefore, in a small enough 
neighborhood of a generic point, a generic field of planes on a manifold of dimension 
three or four is a Goursat structure. Note, however, that starting from dimension 
five the class of Goursat structures is of infinite codimension within the class of all 
rank two distributions. Indeed, for a generic field of planes D on K", for n > 5, we 
have dimT)^ 1 \p) = 3 and dimD^fjo) = 5, for any point p on an open and dense 
subset of R n . 

We give now a mechanical example. For n > 0, the distribution spanned by the 
following pair of vector fields: 

n-l 

J- , cos(#o)7ro^ + sin(^ )vr ^ + E sin(0 i+1 )7r i+1 ^ 

i=0 

where 7Tj = YYj=i+i cos (^j) anc ^ = 1, is a Goursat structure on R 2 x (S l ) n+1 
equipped with coordinates (£i,£2,#0) • • • ,@n)- This distribution is the kinematical 



model for the "nonholonomic manipulator" of S0rdalen, Nakamura, and Chung |63 
Another example is the n-trailer system (see Section |3|), which will play a funda- 
mental role in this article. 

1.3 Goursat Normal Form 

The concepts of derived flag and Goursat structure were introduced, using the dual 
language of Pfaffian systems, by E. von Weber |7(| in order to characterize the 



class of Pfaffian systems that can be converted into a particular normal form, also 
introduced by him, which is now known as Goursat normal form (see (|3]) below; see 
also Appendix |A|) . 

Although it is not clear who was the first to prove that Goursat structures can 
be converted locally into Goursat normal form, at least on an open and dense subset 
70], Theorem VI] (compare @, 0, |2|, g|, [0, and J7|); it is clear that the 



foundations of this result were set by Engel and Weber (see 0, [|T(J, |70| , and 
the references given there; see also Appendix |A|). The importance of this result 
was, however, fully understood only later, by E. Cartan when he solved a long 
standing problem of that time: the characterization of explicitly integrable Monge 
equations Q (see also §, @, [0, H, and III)- 



11 



Theorem 1.2 (Weber-Cartan-Goursat) Let V be a Goursat structure defined 
on a manifold M of dimension n > 3. There exists an open and dense subset 
U C M such that, for any point p in U , the distribution T> is locally equivalent at p 
to the distribution spanned by the following pair of vector fields: 



8x n 



X r. 



dx„ 



-i +Xn - 1 dx n ^ 2 ^ \~ X3 dx 2 + 8X1 



(3) 



considered on a small enough neighborhood V cM, n of zero. Moreover, ifn = 3 or 4 
then U = M. 

In control theory, the normal form ([J) is usually called chained form. In order to 
keep the classical name, we will call it Goursat normal form. An elegant charac- 
terization, using the growth vector, of the above mentioned open and dense set U 
was obtained by Murray |51] (see Theorem |5.2| ). Observe that, in most of the above 
mentioned references, Goursat structures are not defined by distributions but by 
their duals, that is by Pfaffian systems. Note also that many other names have 
been given to Goursat structures: "systeme vom Charakter eins und vom Rang 
zwei" |70[| , "systemes de classe zero" jLlj], "systemes speciaux" f22|, "systemes en 
drapeaux" , "systems of Goursat type" [0] , and "systems that satisfy the Gour- 



sat condition" [49]. 



2 Kumpera- Ruiz's Theorem 
2.1 Kumpera- Ruiz Normal Forms 

If at a given point a Goursat structure can be converted into Goursat normal form 
then this point is called regular] otherwise, it is called singular. The first who 
observed the existence of singular points were Giaro, Kumpera and Ruiz ||20|| . This 



initial observation has led Kumpera and Ruiz to write their pioneering paper pz 
where they introduced the normal forms that we will consider in this section. 

We start with the Pfaff-Darboux and Engel normal forms, given respectively 
on M 3 , equipped with coordinates (xi,X2,xs), and R 4 , equipped with coordinates 
(xi, X2, x 3 , X4), by the pairs of vector fields k 3 = (k?,/^) an d = (k 4 ,^), defined 
by 

K 3 _ 8 
1 dxz 

^2 ^3 q X2 Q X1 
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and 

,A = d 

1 8x4 

Loosely speaking, we can write 

K l ~ dx 4 

In order to make this precise we will adopt the following natural notation. Consider 
a vector field 



on M™ -1 equipped with coordinates (x±, . . . ,x n -i). We can lift / n_1 to a vector 
field, denoted also by / n_1 , on IR n equipped with coordinates (xi, . . . ,x n -i,x n ) by 
taking 

p-i = v-»-ip-i( _ \_a_ + n._a_ 

J 2-<i=l Ji l x l' • • • i x n-l) dx . -T U 9a . n . 

That is, we lift / n_1 by translating it along the avdirection. 

Notation 2.1 From now on, in any expression of the form k% = a(x)K^ 1 + 
(3(x)K2~ l ) the vector fields k™ -1 and k^~ x should be considered as the above de- 
fined lifts of k^ 1 and «2~\ respectively. 

Let K a ~ l = (ft™ -1 , ^2~ l ) denote a pair of vector fields on W 1 ^ 1 . A regular prolon- 
gation, with parameter c n , of n n ~ l , denoted by K n = R Cn (K n ~ 1 ), is a pair of vector 
fields K n = (k™, /eg) defined on R" by 

n _ d 

1 fen (4) 

where c„ belongs to R. The singular prolongation of K n_1 , denoted by /t n = S'(/t™ _1 ), 
is the pair of vector fields K n = defined on M n by 

n _ 

1 dXn , ! (5) 

Definition 2.2 A pair of vector fields n n on M n , for n > 3, is called a Kumpera-Ruiz 
normal form if n n = o- n -3 ° • • • ° o"i(^ 3 ), where each <jj, for 1 < i < n — 3, equals 
either S or R c ., for some real constants q. 
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We will also call a Kumpera-Ruiz normal form the restriction of a Kumpera-Ruiz 
normal form to any open subset of ¥L n . Most often, the coordinates x%, ... ,x n will 
be the elements of a coordinate chart x : M — > M. n , defined in a neighborhood of a 
given point p in M. If we have x(p) =0 then we will say that the Kumpera-Ruiz 
normal form is centered at p. For example, the two models considered in |20] : 



O rp 9 I ~, 9 I ~, 9 I 

9xs ' 5 9^4 * <9x3 *J 9x2 9xi 

9 9 _i_ ™ f „ 9 _i_ ™ 9 I 9 

9x 5 ' 9x 4 ^ X5 \ XA dx 3 ^ X3 9x 2 ^ 9xi 

defined by Rq(k 4 ) and 5"(k 4 ), respectively, are Kumpera-Ruiz normal forms on M 5 , 
equipped with coordinates (x%, . . . , 25), centered at zero. 

2.2 Kumpera-Ruiz 's Theorem 

The following theorem of Kumpera and Ruiz shows clearly the importance of their 
normal forms. We will give a proof of this theorem at the end of this Section because 
many of our results are based on a construction that also appears in our proof. 
Moreover, we would like to emphasize two features of our proof. Firstly, it is quite 
close to the original ideas of E. von Weber. Indeed, though we use distributions 
instead of Pfaffian systems, the two proofs share the same fundamental Lemma 
(compare |7J], Theorem V] and Proposition |275|; see also Appendix |AT). Secondly, 
it is to our knowledge the only one that does not mix the language of vector fields 
and differential forms (everything is done in terms of vector fields). For alternative 
proofs we refer the reader to JT2J and to the original work of Kumpera and Ruiz 



Theorem 2.3 (Kumpera-Ruiz) Any Goursat structure on a manifold M of di- 
mension n is locally equivalent, at any point p in M, to a distribution spanned by 
a Kumpera-Ruiz normal form centered at p and defined on a suitably chosen neigh- 
borhood [/CM" of zero. 

This theorem is the cornerstone to understand the geometry of Goursat struc- 
tures. On the one hand, this result implies that locally, even at singular points, 
Goursat structures do not have functional invariants; this property makes them 
precious but rare and distinguishes them from generic rank two distributions on 
n-manifolds, which do have local functional invariants when n > 5 (see JTO , ]2"7 |, 



68], [76]). But on the other hand, the real constants that appear in Kumpera-Ruiz 



normal forms are unavoidable; this fact has been observed only recently and implies 
that Goursat structures do have real invariants (see JIB] , p9j, and Section |7|). 
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Though our definition of Kumpera-Ruiz normal forms was inductive, it is possible 
to give the following equivalent explicit definition (^|), which will also be used in 
the paper. Observe that in the normal form (|6]), we use a double indexation x*- of 
coordinates, for < % < m+1, where the integer m gives the number of singularities 
of the normal form, that is the number of singular prolongations (provided that o\ 
is regular, which can always be assumed without lose of generality). 

Corollary 2.4 Any Goursat structure defined on a manifold M of dimension n > 4 
is locally equivalent, at any point p in M , to a distribution spanned in a small 
neighborhood of zero by a pair of vector fields that has the following form: 

m i—1 ki — l 

(aSj , E(IK)(E (4 + 4)^7 + ^) )> ( 6 ) 

1 i=o j=o j=l J+ 1 u i 

where the coordinates Xp for < % < m + 1 and 1 < j < hi, are centered at p; the 
integer m is such that < m < n — 4; and the integers hi, for < % < m + 1, satisfy 
k > 1, . . . , A; m _i > 1, k m > 3, k m+1 = 1 and ^ = n - The constants c*, for 

1 < j < A;j — 1, are real constants. 



2.3 Low Dimensional Examples 

Let us recall the complete local classification of Goursat structures on manifolds of 
dimension five and six, obtained by Giaro, Kumpera and Ruiz (see and |32[|). 

(i) Any Goursat structure on a five-manifold is locally equivalent to one of the 
two following Kumpera-Ruiz normal forms 



( a§7 ' X5 afl + ^afi + X3 afl + air ) ^ 

( at ' at + Xs ( X4 A + x ^ + ) ' ( 8 ) 

which are not locally equivalent at zero. 

(ii) Any Goursat structure on a six-manifold is locally equivalent to one of the 
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five following Kumpera-Ruiz normal forms 



9 

8x 6 

a 

8x6 

8 

8x 6 

a 

9x6 

d 
8x6 



X& 9x 5 + Xs 9x 4 + Xi 8x 3 + Xs dx 2 + 9xi 



8x4 

c 

8x4 



8x3 

c 
8x3 

8 
8x3 



8 _J_ I J. 8 I ~, 8 I rp 8 I 9 

9xr 6 15 Qj.. ~r 4 q x ~t~ 3 g ~ q 



Xe 9x 5 + 9x4 + X5 \ Xi 9x 3 ~^ Xs 9^2 ~^ 9^ 
(*6 + 1) ^ + all + *5 + + £ 



9 
9x5 



(9) 
(10) 

(11) 
(12) 

(13) 



which are pairwise locally non-equivalent at zero. Observe that these two results are 
not implied by Theorem |2.3| . Indeed, they show that in dimension five and six the 
constants that appear in Kumpera-Ruiz's Theorem can be "normalized" to either 
or 1. 

For n = 7, 8 and 9 the complete classification is more delicate, but there is still 
a finite number of models (see Rl2"| , |T9|| , [0, and p9|). For n > 10, the number of 
local models is infinite (see [4~9"|, and Section |7]) and the complete classification 
remains an open problem (see recent results in ||50||). 



2.4 Proof of Kumpera-Ruiz's Theorem 

Our proof of Theorem pT3| will be based on the following classical result, which was 



originally formulated in the dual language of Pfaffian systems |70| , Theorem V] (see 
also Q, @, |32|], and Appendix 0). 



Proposition 2.5 (E. von Weber) Let T> be a Goursat structure on a manifold 
M of dimension n > 4. In a small enough neighborhood of any point p in M , the 
distribution T> is equivalent to a distribution spanned on M. n by a pair of vector fields 
that has the following form: 



( 8 8 y n ' VnCl 



n— 1 I f-n— 1 



(14) 



where Q_ 1 and Q 1 are the lifts of a pair of vector fields that span a Goursat 



structure on 



DTI— 1 



and the coordinates yi, ■ ■ ■ ,y n are centered at p. 



Proof of Proposition |2.5| It is well known (see e.g 



and |75|] ) that any 

Goursat structure Dona manifold of dimension n > 4 admits a canonical line field 
C C V> uniquely defined by [C, T>^'\ C T>^\ Observe that in the preliminary normal 
form flT4j ) of Proposition ^7S] this line field is given by C = (^-). 
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It is clear that, applying around p the flow-box theorem to a vector field that 
spans £, we can chose local coordinates (zx, . . . , z n ), centered at p, such that V is 
locally equivalent to a distribution spanned on W 1 by a pair of vector fields that has 
the following form: 




where C = (^f-)- Since dimZ>M(p) = 3 there exists an integer i such that ff^G ) 
0. We can assume that i = n — 1 and, moreover, that a„_i(0) = 0. Otherwise, 
replace the coordinate z n _i by z n _i — Zia: n _i(0). Now, if we define y n = a n _i(z) 
and yi = Zi, for 1 < i < n — 1, we get that V is locally equivalent to a distribution 
spanned on R n by a pair of vector fields that has the following form: 

( air ' ^a^hr + ^1=2 ^(y)wi + Wi ) ■ 

But the inclusion [£,£>W] C £> (1) clearly implies ^ = for 2 < i < n - 2. That 
is f3i(y) = a i (y n _ 1 )y n + &i(y n _i), where ^ n _ x = (y 1? . . . ,y n _i). Define 

cr^^tr + E&^-i)!- and cr 2 = E::My n -i)i- + i- 1 - 

We conclude that V is equivalent to ( ^ , ^C^+C:?" 1 ), where both C™ -1 and C^" 1 
are lifts (see Notation [Hp of vector fields defined on W 1 ' 1 . Put T = (Ci™ -1 , C 2 ™~ 1 )- 
Clearly dimV^ +1 ^ = dimjF^) + l, for < i < n—3. It follows that the distribution T 
is a Goursat structure on R n_1 . □ 



Proof of Theorem |2.3| On three-manifolds, Theorem ^73] is obviously true. Indeed, 
it is the solution of the Pfaff problem (see [15] and [18]; see also || for a modern 
approach). We will proceed by induction on n > 4, showing that if any Goursat 
structure on an (n — l)-manifold can be converted locally into a Kumpera-Ruiz 
normal form then the same is true for any Goursat structure on a manifold of 
dimension n. 

Let D be a Goursat structure on a manifold M of dimension n > 4 and let p be 
an arbitrary point in M. It follows from Proposition |2.5| that T> is equivalent, in a 
small enough neighborhood of p, to a distribution spanned on M. n by a pair of vector 
fields (£", that has the following form: 



si 

/-n 
S2 



a 

VnCl 



n— 1 1 f-n— 1 
T S2 
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In the rest of the proof we will assume that T> = Q )• Note that the ^/-coordinates 
are centered at zero. 

Our aim is to build a local change of coordinates 

Oi,... ,x n ) = ( j) n (y 1 ,... ,y n ), 

a Kumpera-Ruiz normal form (k^k^) on MJ 1 , and three smooth functions /i n , v n , 
and 7] n such that 

<t>m) = (v n orK (15) 

where ip n = (0™)™ 1 denotes the inverse of the local diffeomorphism <p n and both 
A i n(0) 7^ and ^ n (0) 7^ 0. Moreover, we will impose the x-coordinates to be centered 
at zero. That is 0" (O) = 0. Observe that the triangular form in (|T5| ) appears because 
both Ci and span the canonical line fields of the distributions spanned by Q) 
and ( 2), respectively. 

By Proposition |2.5| , the distribution spanned by (£™, C^) is defined by the lifts of 
a pair of vector fields (Ci" _ \ C^ -1 ) that span a Goursat structure on IR n_1 . Since the 
Theorem is assumed to be true on M n_1 , the distribution spanned by (CT~\ C2 _1 ) i s 
locally equivalent to a Kumpera-Ruiz normal form ( ) defined on IR n 1 and 

centered at zero. It follows that there exists a local diffeomorphism (xi, . . . , x n -i) = 
. . . ,y n -i) an d four smooth functions z/ n _ 1; A n _i, r] n _i, and [i n -\ such that: 



(16) 



where tp n 1 = (0™ *) 1 denotes the inverse of the local diffeomorphism <p n 1 and 
(^n-lyUn-l - A„_ir? n _i)(0) 7^ 0. 

Let (p n = (0 n_1 , 4> n ) T be a diffeomorphism of M. n such that depends on the 
first n — 1 coordinates only. Moreover, let / be a vector field on R n of the form 
/ = a/™ -1 + f n , where a is a smooth function on MP, the vector field / n_1 is the 
lift of a vector field on R n_1 (see Notation p.l| ), and the only non-zero component 
of f n is the last one. A direct computation shows that: 

€(/) = (a o rW^iP- 1 ) + ((L/0„) o ^) (17) 

Note that the vector field 0™ _1 (/ n_1 ) is lifted along the x„-coordinate, which is 
defined by <p n . 
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Regular case: If /i n -i(0) 7^ then we can complete <p n 1 to a zero-preserving local 
diffeomorphism of M. n by taking <fr n = (<p n ~ 1 , <f) n ) T , where 

, / v "n-lVn + Vn-l Vn-l(0) 
<Pn(y) = T T 77TT- 

In this case, we define c n = (^„_i//U n _ 1 )(0) and 

v n = LqL(/) n , 7] n = hQ(j) n , and fi n = \ n -iy n + ju n -i- 

Observe that z/„(0) = L<^0„(O) = (z/ n _i/i„,_i — A n _i?7 n _i) (0) 7^ and that /x n (0) = 
/i n _i(0) 7^ 0. Thus the right hand side of (|15|) defines a locally invertible trans- 
formation. Moreover, the Kumpera-Ruiz normal form (/«", k^) is defined to be the 
regular prolongation, with parameter c n , of (ft™ -1 , «2~ )■ 

Let us check that, in this case, relation fll5| ) holds. Together, relations ([16|) 
and flT7D give: 

c(c 2 n ) = (y„ o r^r^cr 1 ) + c^cr x ) + ((l^) ° r) £ 

= {{Vn-lVn + Vn-l) ° V*") ^ + {{K-lVn + l*n-l) ° 4" 1 + ° 

= ((v- iy „ + Mn-o o r) (( z:ZZ:~-\ ° <~ a + + fo» ° 

= o tfj n ) ((x n + Cn) 4" 1 + KT 1 ) + (Vn o VKl 

Moreover, we have 

c(cd = ((L cr 0„) of)i-, [v n o ^k+ 3 . 

It follows that, in the regular case, relation ([TJ|) holds. 



Singular case: If /x n _i(0) = then we can complete n 1 to a zero-preserving local 
diffeomorphism of M™ by taking n = (0 n_1 , 4> n ) T , where 

<Pn(y) = ; • 

Pn-lVn + Vn-l 

Observe that /U n _i(0) = implies 0„(O) = 0. Additionally, since /i„_i(0) = and 
[y n -\^n-\ - A n _ 1 ry n _i)(0) 7^ 0, we have A n _i(0) 7^ and r/ n „i(0) 7^ 0. In this case, 
we define 

v n = L C n0 n , r] n = L C n(f) n , and /i n = v n _ x y n + r/ n _i. 
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Observe that i/ n (0) = L C n0 n (O) = (i/ n -i/V-i - A ra _i77„_i)(0) 7^ and that /i n (0) = 
77^—1 (0) 7^ 0. Thus the right hand side of fll5|) defines a locally invertible trans- 
formation. Moreover, the Kumpera-Ruiz normal form (k™, k^) is defined to be the 
singular prolongation of («" , ^2 _1 )- 

Let us check that, again, relation (15]) holds. Together, relations fll6l) and fll? ) 



give: 

= ((^-iy n + »7n-l) + ((A n -U/„ + ^-l) ° ^2 ^ + fan ° f)«J 

= (Mn O V") (AC?" 1 + X n «5" X ) + (t?„ o 

Like in the previous case, we have 



It follows that relation ([15]) holds in both cases. □ 



3 The N- Trailer System 



3.1 Definition of the N- Trailer System 

The kinematical model for a unicycle-like mobile robot towing n trailers such that 
the tow hook of each trailer is located at the center of its unique axle is usually called, 
in control theory, the n-trailer system (see the books |36| and |38| ; the papers [17| , 



25 , 29 



. |66| , and ; and the references given there) . Figures 
representing this system are given in Appendix 0. For simplicity, we will assume 
that the distances between any two consecutive trailers are equal. 



Definition 3.1 The n-trailer system is the distribution defined on R 2 x (S l ) n+1 , 
for n > 0, by the following pair of vector fields: 



^- , vr cos{V ) JL + vr sin(0 o ) ^ + E fi+i sin(^+i - 0*) J- 



n-l 

E< 

i=0 



;is) 



where 7r,- 



and 7r n = 1. 
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In the above definition, the functions £1, £2, 0o>---> #n axe coordinates on the 
manifold M 2 x (S ,1 )™ +1 . The coordinates £1 and £2 represent the position of the 
last trailer, while the coordinates 60, . . . ,6 n represent, starting from the last trailer, 
the angles between each trailer's axle and the £i-axis. It is easy to check that the 
n-trailer system is a Goursat structure (see e.g. [j35H ). 

We give now an equivalent inductive definition of the n-trailer. This defini- 
tion already appears in |2{J and reminds the one given in the previous section for 
Kumpera-Ruiz normal forms. To start with, consider the pair of vector fields (r°, t®) 
on M 2 x 5* 1 that describe the kinematics of the unicycle-like mobile robot towing no 
trailers: 



^0 _ d 
r i ae 



r 2 = cos(^o)4 + sin(^o)^ 



d_ 
2 ' 



The n-trailer system can be defined by applying successively a sequence of prolon- 
gations to this mobile robot. In order to do this, suppose that a pair of vector fields 
r n_1 = (t" -1 ,^ 1-1 ) on M 2 x (S 1 ) 71 , associated to the mobile robot towing n — 1 
trailers, has been defined. The pair of vector fields r n = (r™, r^) on M? x (S' 1 ) n+1 
corresponding to the n-trailer system is then given by 



T n = -2- 

1 ae„ 



Observe that this definition should be understood in the sense of Notation [2.1| , and 
that the pair of vector fields (rf , r^) coincides with that of Definition p.lj Mechan- 
ically, to prolongate the n-trailer means to add one more trailer to the system. 



3.2 Local Conversion of the N-Trailer System into a Kum- 
pera-Ruiz Normal Form 

Since the n-trailer is a Goursat structure, it follows directly from Kumpera-Ruiz's 
theorem that, in a small enough neighborhood of any point of its configuration 
space, in particular at any singular configuration, the n-trailer can be converted into 
Kumpera-Ruiz normal form. In this subsection we describe this conversion explicitly. 
For regular configurations, our result gives the transformations proposed in |5Dfl 
and [|57[] ; for singular configurations, our result gives a new kind of transformations. 



Denote by ( the coordinates of M 2 x (S ,1 ) n+1 , that is 

C = (Ci> — >Cn+3) = (€l,£2,@0, —,0 n ). 



21 



Fix a point p of R 2 x (S 1 ) 11 ^ given in (-coordinates by £(p) = C p = (£1, £2, C 

In order to convert, locally at p, the n-trailer into a Kumpera-Ruiz normal form we 

look for a local change of coordinates 

Ol, . . . , X n+3 ) = 0"(£i, £ 2 , 00, ■ ■ ■ , On), 

a Kumpera-Ruiz normal form (fti +s , k^" 1 " 3 ) on M n+3 , and three smooth functions u n , 
r) n , and /i„ such that 

W(T?) = (l)„ofK +3 + ( fr of) fi f, 

where ^ n = (0™) -1 denotes the inverse of the local diffeomorphism 0™ and both 
u n(C p ) an d yUn(C p ) 7^ 0- Observe that we do not demand the x-coordinates to 
be centered at p, and thus the point x{p) = (0 n o ()(p) will be, in general, different 
from zero. 

We construct inductively here a change of coordinates n = (0i, . . . , n+3 ) T and 
three smooth functions v n , r) n , and p, n that satisfy (|TP|). We start with n — 0. If 
6*q 7^ ±7r/2mod27r then define 0i = £i, 02 = 62, an d 03 = tan(# )- Moreover, take 

= cos(6'o), Vq = sec 2 (6'o), and r/ = 0. If 6>q = ±7r/2mod27r then define 0i = £ 2 > 
02 = £1, and 03 = cot(6 l o). Moreover, take /io = sin(6 l o), vq = — csc 2 (6*o), and r] = 0. 
Denote = sin(#j — and Cj = cos(#j — for < % < n. 

Now, consider the sequence of smooth functions defined locally around the 
point C{p)i f° r 1 < 2 < n, by either 



l^i Q/^i— 1 

v i — L T i0j + 3 

if #f — 0f_ x 7^ ±7r/2mod27r (regular case) or by 



CilH-\ 

(20) 



■>i+3 



//i = SiVi-i + Q?7j_i 
= L T i 0j + 3 

»7» = L T »0 i+3 , 



(21) 
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if Q\ — 6 l f_ 1 = ±7r/2mod27r (singular case). It is easy to prove that, for < i < 
n, the transformations defined by (x±, . . . ,^+3) = £2, @o, ■ ■ ■ are smooth 

changes of coordinates around Pi and that, moreover, we have both ^(C Pi ) 7^ and 
/ii(C p< ) 7^ 0, where Pi denotes the projection of p on M. 2 x (5 a ) l+1 , the product of IR 2 
and the first i + 1 copies of S\ that is C Pl = (£?, & 0g, • • • , 0?). 

Proposition 3.2 For n > 0, the local diffeomorphism <p n and the smooth functions 
v n , r] n , and fi n satisfy fllTf), and i/ras convert locally the n-trailer system into a 
Kumpera-Ruiz normal form. 



Proof of Proposition |3.2| We will prove that the relation (O) holds for n > 



by induction on the number n of trailers. Relation ( |19D is clearly true for n — 0. 
Assume that it holds for n — 1 trailers, that is 

0r 1 (rr 1 ) = K-io^- 1 x 2 

The inductive definition of the n-trailer gives 

a 

_n— 1 1 __„//) /) \_n— 1 



r 2 " = sin(#„ - O^rr 1 + cos(0„ - 6^)1% 
Recall (see the Proof of Theorem |2l|) that for a diffeomorphism n = (0 n_1 , <pn+^) T 



of M n+3 , such that n_1 depends on the first n + 2 coordinates only, and for a vector 
field / on M n+3 of the form / = af n ~ l + /n+3, where a is a smooth function on M n+3 , 
the vector field f n ~ l is the lift of a vector field on M n+2 (see Notation |2.1| ), and the 
only non-zero component of f n+ % is the last one, we have 

€(f) = (« o r)<p:~\r- 1 ) + ((L/0„+3) o r) ^ (22) 

Note that the vector field 0*~ 1 (/ ?1 ~ 1 ) is lifted along the x n+3 -coordinate, which is 
defined by 4> n+3 . 

In the regular case, that is if 9f — Q\_ x ^ ±7r/2mod27r, we take a regular pro- 
longation and, using relations fl20|) and (|22|) , we obtain: 

w (r«) = ( Sn .o^)0r i (rr i ) + (c n o^)0r i (rr i ) + ((L T n0 n+3 )o^) " 



*K'2J — K°n^<y JY* V'l )^\^n^<Y JY* V'2 i T " f 

= (( Sn l/ n _! + Ctfn-l) o V H ) < +2 + ((Cn/in-i) o tfj n ) K% +2 + ( Vn o ^K +3 

= (/in O (*„ +3 < +2 + 4 +2 ) + for, ° ^"K +3 

= feof)< +3 + (p n of)4 1+3 
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In the singular case, that is if #f — 8f_ l = ±7r/2 mod27r, we take the singular 
prolongation and, using relations ( pi]) and (22), we obtain: 

+ (r /n o^)<+ 3 



+.3 



Moreover, in both cases, we have: 



((L rr n+3 ) o ^ = (i/ n o rx 3 - 



It follows that relation ( |i9|) holds for any n > 0. 



□ 



3.3 Local Conversion of an Arbitrary Goursat Structure 
into the N- Trailer System 



Reversing the construction given in the Proof of Proposition [O] leads to the follow- 
ing surprising result (already announced in [54) and proved in fl55|), which states 



that the n-trailer system is a universal local model for all Goursat structures. This 
theorem will play a fundamental role in this article. Indeed, it will allow us to gen- 
eralize local results known for the n-trailer, like the formula for the growth vector 



obtained by Jean |29|, to all Goursat structures. 



Theorem 3.3 Any Goursat structure on a manifold M of dimension n + 3 is locally 
equivalent, at any point q in M, to the n-trailer considered around a suitably chosen 
point p of its configuration space M. 2 x (S 1 ) n+1 . 

Proof of Theorem |3.3| By Theorem |2.3| , our Goursat structure is, in a small 
enough neighborhood of any point q in M, equivalent to a Kumpera-Ruiz nor- 
mal form K n+3 . Denote by y = (yi, . . . ,y n +3) the coordinates of K n+3 and put 

ivl--- ,2/n+s) =y(?)- 

Recall that, by definition, the pair of vector fields n n+3 is given by a sequence of 
prolongations n % = 0V3 o • • - o cr^/t 3 ), where <7j belongs to {R Cj , S}, for 1 < j < i — 3 
and 3 < i < n + 3. We call a coordinate yi such that k % = S^k* -1 ) a singular 
coordinate, and a coordinate yi such that k % = R c (k 1 ~ 1 ) a regular coordinate. It 
follows from the proof of Theorem |2.3| (see Section |2|) that for all singular coordinates 
we have y\ = 0; but for regular coordinates, the constants y\ can be arbitrary real 
numbers. 
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To prove Theorem |3.3| , we will define a point p of M 2 x (S' 1 )" +1 whose coordinates 
C(p) = C p = (£1 , £f, ^o) • • • > d n) satisfy (x o £)(p) = y(q), where x and C denote the 
coordinates used in the Proof of Proposition p^2| . First, put the axle of the last 
trailer at (2/15 2/2) 5 that * s (£15 £2) = iViiVV)) an d take 9$ = arctan(yf). Compute 
X3 = tan(#o), Ho = cos(6 l o), Vq = sec 2 (9 ), and ^ = 0. Then, take for % — 1 up to n, 
the following values for the angles 8f mod 2ir. If the coordinate yi + s is singular then 
put 6? = Q\_ x + 7r/2 and compute the coordinate Xj +3 and the smooth functions /ij, 
z/j, and r/i using (|2~lD . If y i+3 is regular then put 



and compute the coordinate 0^+3 and the smooth functions 24, and 773 using (p0|) . 
The result of this construction is that (x o C)(p) = £/(<?)• By Proposition |3^2| , the 
coordinates io( convert the n-trailer into a Kumpera-Ruiz normal form. By the 
above defined construction, this normal form has the same singular coordinates 
as /t n+3 and is defined around the same point of R n+3 (if we translate the regular 
coordinates in order to center them then those two Kumpera-Ruiz normal forms 
have the same constants in the regular prolongations). Hence, the diffeomorphism 
o x^ 1 o y gives the claimed equivalence. □ 



4 Singularity Type 

4.1 Characteristic Distributions 

A characteristic vector field of a distribution T> is a vector field / that belongs to 
T> and satisfies [/, T>\ C T>. The characteristic distribution of a distribution T> is 
the module, over the ring of smooth functions, generated by all its characteristic 
vector fields. A characteristic distribution need not be of constant rank; but it 
follows directly from the Jacobi identity that a characteristic distribution is always 
involutive. The main ingredient in the definition of the singularity type will be the 
characteristic distributions Cj defined by the following result, which is apparently 



due to Cartan [JTTJ , although he did not state it explicitly in his published works. 
Its proof can be found in |!2| and |43[] (see also 0, and Appendix |A|), were 
slightly stronger versions are proved using the dual language of Pfaffian systems. 

Proposition 4.1 (E. Cartan) Consider a Goursat structure T> defined on a man- 
ifold of dimension n > 4. Each distribution forO <i< n—A, contains a unique 
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involutive subdistribution d C that is characteristic for and has constant 

corank one in 2?^ . 



4.2 Jacquard's Language 

An alphabet is a finite set A whose elements are called letters. A word over A is a 
finite sequence w = (wx, . . . ,Wi), where each u>, belongs to A for 1 < i < I. The 
integer / is called the length of the word w and we denote it by \w\. The empty word 
is the only word of length 0. We denote it by e. By abuse of notation, we will often 
write W\- • - Wi instead of (wx, . . . , w{) and a 1 instead of (a, . . . , a), for any letter a 
repeated I times. Now, call A* the set of all words, including the empty word, over 
the alphabet A. A language over A is a subset L of A*. The concatenation of two 
words v and w over A is the word vw = (vx, ... ,V[, Wx, . . . , u> m ), where I — \v | and 
m — \w\. The concatenation of a language L and a word u> over A is the language 

Lw = {uw : u E L}. 

The shift of a word w = (wi, . . . , uv) is the word (it;)' = (wx, . . . , Wj-i). By defini- 
tion, we take (e)' = e. 



We define now a sequence of languages that will play an important role in this 
paper, since they will describe all possible singularity types of a Goursat structure. 
For a fixed n > 0, consider the alphabet A n = {a , . . . , a n _x} (note that A = 0). 
Following El define, inductively, the Jacquard language J n by J = {e}, Ji = {a }, 



and 

■Al = Jn-l a U Jn-l a l U ^n-2&lfl2 U . . . U JxCLxO-2 ' ' ' &n-l) 

for any integer n > 2. Clearly, the elements of J n are words over A n that all have 
length n. For example, we have J2 = {aoaoj^o^i} and 

J3 = {^O^O^O) ttofllflO) ^oOlfll) CI0CI1CI2}. 

It is easy to check, using an induction argument, that 

card(J n ) = 3 card(J n _i) — card(J„_ 2 ) ) 



for n > 2 (see \2§) 
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4.3 Definition of the Singularity Type 

We start with the definition of a sequence of canonical submanifolds, which will 
lead to the definition of the singularity type. Let V be a Goursat structure on a 
manifold M of dimension n. For < i < n — 5, define the subset C M by 

S® = {peM:V®(p)=C i+1 (p)}, (23) 

where the distribution Ci denotes the characteristic distribution of (see Propo- 

sition (4.1p. For i > n — 4 define = 



Furthermore, starting from S'q 4 -' define, for 1 < j < i, a sequence of subsets 
M D Sf D • • • D Sf x D Sf D • • • D sf by 

S® = { P e sjl : P^&O n T^«! ^ ^(p)}- (24) 

For j > i + 1 define = 0. The above definition is correct since, for any non- 
negative integers % and j, the subset C M is a smooth submanifold of M (we 
consider an empty set as smooth). Indeed, we have the following result, which will 
be proved later, in Subsection [4.5|. 



Proposition 4.2 Let T> be a Goursat structure on a manifold M. 

(i) For any non-negative integers i and j , the subset Sj C M defined by the 
relations (^3]) and is either empty or a smooth submanifold of M that has 
codimension j + I in M. 

(ii) For any non-negative integers i, j and k such that k ^ j we have the following 
relation: Sf +j) fl = 0. 

The fact that a point p belongs to some submanifolds S^, for two given non- 
negative integers % and j, is invariantly related to the Goursat structure at this 
point p. This information, however, is in general redundant. For example, if p 
belongs to Sj we know, by definition, that p belongs also to Sj_ k , for < k < j, 
and that, by Proposition |4~2"1 , it does not belong to Sj l ^\ for k ^ 0. In the following 
definition we propose a way to encode the essential part of this information in a word 
of the Jacquard language (see Corollary [D| at the end of Subsection |4 . 5| ) . 



Definition 4.3 Let T> be a Goursat structure defined in a neighborhood of a fixed 
point p in a manifold of dimension n. The singularity type of T> at p is the word 

5 v {p) =w n _ 4 ---w 

defined as follows: For < i < n — 4, we take Wi = a J+ i if there exists some integer 
j > such that p belongs to we take Wi = ao otherwise. 
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The correctness of this definition is assured by Proposition |4.2| , which states 
that if there exists an integer j > such that p belongs to then this integer is 

unique. For some readers this definition may seem rather abstract at a first glance; 
but it appears quickly, once computed for concrete Goursat structures, that the 
singularity type really encodes essential geometric information. In fact, we will see 
that the growth vector and the abnormal curves of a Goursat structure are given by 
its singularity type (see Sections | and |6[). 

The singularity type should not be confused with the codes used in |L2] and |49|] to 
label Kumpera-Ruiz normal forms. Indeed, the singularity type is, by construction, 
an invariant of a given Goursat structure; while the codes of EI and [03 are not 
invariant: a given Goursat structure can have different codes, depending on how it 
is represented by a Kumpera-Ruiz normal form. 



4.4 Low Dimensional Examples 

For any Goursat structure on a manifold of dimension three or four the singularity 
type is equal, at any point, to e or ao, respectively. That is, the singularity type of a 
contact or an Engel structure does not depend on the point at which the distribution 
is considered. This should be compared with the singularity type of a Goursat 
structure on a five-manifold, which can be either a^ao or a^a\ at a given point p, 
depending on whether or not the Goursat structure can be converted into Goursat 
normal form in a small enough neighborhood of p. Indeed, for the Goursat structure 
spanned by the regular Kumpera-Ruiz normal form (]?]) the canonical submanifold 
Sq is empty, and thus the singularity type equals a a at each point of M 5 ; for the 
Goursat structure spanned by the singular Kumpera-Ruiz normal form (|8|) we have 
^'q '* = {x$ = 0}, and thus the singularity type equals a ai if x$ = 0; and aoa if 
x 5 ^0. 

Let us give one more example, in dimension six. Consider the distribution T> 
spanned by the following Kumpera-Ruiz normal form on R 6 : 




where the constant c$ is either equal to or 1. For both values of Cq, we have 
= and = {x 5 = 0}. For c 6 = 1 the submanifold S± is empty (in a small 
enough neighborhood of zero); for cq = we have = {x§ = Xq = 0} . Therefore, 
the singularity type of V at zero equals a aia if c 6 = 1; and equals a aia 2 if c 6 = 0. 
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4.5 Singularity Type of Kumpera-Ruiz Normal Forms 

Let k u be a Kumpera-Ruiz normal form on W 1 obtained by the inductive procedure 
described in Section Q. We define inductively the word 5 K n of J n _ 3 by 5 K 3 = e, 
5« 4 = clq, and, for n > 5, by the relations 

o~n n = S K n-i ai if K n = S'(k™~ 1 ); 

& = S K n-i aj + i if K n = Ro(k u ~ 1 ) and 5^-1 = 5 K n-2 cij, % > 1; (25) 
^k" = ^k™- 1 a o otherwise; 

where the maps S and i? c denote respectively the singular and regular prolongations 
defined in Section |2|. This definition leads to a characterization of the singularity 
type in the coordinates of the Kumpera-Ruiz normal form (see Corollary [4.6| be- 
low). We start with a Proposition that will allow us to prove the results that were 



announced, without proof, in Subsection 4.3 



Proposition 4.4 Consider the distribution defined on ~R n by a Kumpera-Ruiz nor- 
mal form K n . For any pair of integers i and j such that < j < i and for any pair of 
words wi and u>2 such that w = wiaia 2 ■ ■ ■ aj + iW2 belongs to J„_3 and |u> 2 | = % — j , 
we have 5 K n = w if and only if zero belongs to . Moreover, if a submanifold 



Sj contains zero then, in the coordinates (x\, . . . , x n ) of the Kumpera-Ruiz normal 
form k u , we have 

Sj ^ = {Xn~i = 0, . . . , X n -i+j = 0}, (26) 

where < i < n — 5 and < j < i. 

The following Lemma is a direct consequence of the definition of Kumpera-Ruiz 
normal forms given in Section (see Definition |2.2|); its proof is straightforward. 



Note that below all vector fields k% 1 should be considered as vector fields on R n , 
obtained from k^ -1 by i successive lifts, as defined by Notation |27i] (see Section §). 

Lemma 4.5 Let D be a Goursat structure on M. n spanned by a Kumpera-Ruiz 
normal form n n = (k", k^). Suppose that K n = <7„_3 o ••• o <7i(/t 3 ) and denote 
K n ~ l = a n -3-i o ■ ■ ■ o cri(/t 3 ). The derived flag ofV is given by 

P« = (^,...,^-, K r), for0<z<n-3. (27) 

The characteristic distribution C{ C ofT>( t+1 ^ is given by 

ft = (&,..., sfc), forO< l <n-A. (28) 
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Moreover, if S K n = W\a\a2 • • • a^awi, where \vj%\ — i — j and a G {ao, a\, dj+i}, then 
we have 



K 



- Sfc=l 3; n-i+fc &n _ i+fc _ 1 + 9Xn _ i _ 1 +Xn-iK-2 % > ( 2 9) 

for < Z < j — 1 . 



Proof of Proposition |4.4| For any fixed integer i such that < i < n — 5, we will 



prove the Proposition by induction on the integer j, for < j < i. We start with 
j = 0. Assume that 5 K n = W\aw2, where ju^l = i and a G A n _i. If a a\ then the 
definition of 5 K n implies 



K 2 ' — (%n—i 4" C n —i) g x _._ ~\~ ^2 ) 

and it follows from relation ( p7[ ) that 

_ f_8_ 9 fV . _i_ r 1 d I n-i-l\ 

Together with relation (^) this expression implies that, for any point p G M n , we 
have V^(p) ^ C i+ i(p), which implies that Sq^ is empty In particular ^ Sg . 
Otherwise a = a 1; and then the definition of 5 K n implies 

n-i _ d I n-i-1. 



and it follows, again from relation (|27|) , that 

Hence, for any point p G M n , we have "D^(j>) = Cj + i(p) if and only if x n _i(p) = 0. 
In particular G Sq . It follows that Proposition [O] is true when j = 0. 
Now, assume that Proposition JO is true up to the integer j — 1 and that 



5 K n = wiaia2 ■ ■ ■ ajaw2i 
where |ty 2 | = i — j and a G {a 1; a , Since 

we have 
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Observe that relation (p^) gives 

^ = (^•••^>«r +i ) ) (30) 

Moreover, it follows from relation (^), taken for I = j — 1, that 

K 2 +J = Ylk=l X n-i+k dx n ^ i+k ^ + dx^-i-i + X n-i^2 * ^ (31) 

If a = ax then the definition of S K n implies 

n—i+j _ 9 I n—i+j— I 



dx 



n — i+j- 



- + X n _j_|_jK 2 



Therefore the vector field k 2 * +,? that appears in the relation (|30J) is, using (31) 
given by 



n-i+j _ d I /p'-'r g 1 5 L /y .^R-i-lN 

ft 2 - aB„_ 4+i _i +J '"-«+jlZ,t=i J -n- 1 +fca I „_, +t _ 1 



Hence, for any p G Sj_i, we have ^(p) fl TpS 1 ]^ = C^j(p), which implies that 
Sj is empty. In particular ^ Sj l \ If a = a then the definition of 5 K n implies 

n—i+j / | \ Q , n—i+j — I 

K 2 — \ X n-i+j -T C) Q^—^—l + K 2 > 

where c ^ 0. Therefore the vector field k^~ 1+:> that appears in the relation ( p0[) is, 
using fl3ll) , given by 



n-i+j _ r \ A 9 i yJ'^ ■ , 9 i 8 I -r 



Thus for any p G we have V^~^{jp) fl TpS 1 ]^ = Ci_j(p), which implies that 
is empty (at least in small enough neighborhood of zero). In particular ^ Sj . 
Finally, if we have CL — &j+l 

then the definition of 5 K n implies 



jn—i+j d i n—i+j — 1 

i+3- 



O I ' 

x n - i+j dXn — — + K 2 



Therefore the vector field k 2 j that appears in the relation (|30D is, using (^T] 
given by 

n-i+j _ sr^j d , .9 , T „n-i-l 

K 2 ~ 2^k=l X n-i+k 8x n -i +k -i + &r„-i-i -t-J-n-*«2 

Thus for any p G we have V^~^(p) fl TpS^ 7^ Ci-j(p) if and only if 

x n _i(p) = • ■ ■ = x n _ i+j (p) = 0. 
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In particular, we have G Sf\ which ends the proof. □ 



Proof of Proposition 4.2 Consider a Goursat structure T> defined on a smooth 



manifold M of dimension n. First Item: Item (i) of Proposition [4.2| follows directly 



from Proposition W\ which states that, in the coordinates of a Kumpera-Ruiz 



normal form, the restriction of each set to a small enough neighborhood of zero 



is either empty or smooth. Indeed, by Theorem |2.3|, the Goursat structure T> is 
locally equivalent, at any point p in M, to a Kumpera-Ruiz normal form centered 
at p; and hence the restriction of each set to a small enough neighborhood of 
any point p in M is either empty or smooth. This obviously implies that the whole 
set Sj 1 ^ is either empty or smooth. 

Second Item: We will prove Item (ii) by contradiction. Let Dbea Goursat structure 
such that at a point p we have p G <Sj* D S% for a given triple of non-negative 
integers i, j, and k, such that k < j. Take a Kumpera-Ruiz normal form k 71 , 
centered at p, that is locally equivalent to T> at p. In the coordinates of K n , we 
haveOG^ +j) n4 i+A;) . Let w = 5 K n be the word uniquely attached to n n by (p5|), 



and denote the letters of w by w = u>„_4 ■ ■ - w . By Proposition |44], we have both 
Wi + i = a.j-i+1, for < / < j, and Wi + i = at-i+i, for < / < k. In particular, we 
have Wi = aj + i and Wi = a k+ i, which is impossible since the inequality k < j implies 
that a k+ i ^ □ 

The three Corollaries listed below follow directly from Proposition and from the 
definition of the singularity type. 

Corollary 4.6 Let T> be the Goursat structure spanned on W 1 by a Kumpera-Ruiz 
normal form k" . We have 

5t>(0) = S K n, 

that is the singularity type at zero of K n is given by S K n . 



Corollary 4.7 Let V be a Goursat structure defined in a neighborhood of a fixed 
point p of a manifold of dimension n. For any integers such that < j < i, the 
point p belongs to Sj if and only if the singularity type of V at p is of the form 
Sx>(p) = W\a\a2 ■ ■ ■ aj + iW2, with |w2 1 = % — j ■ 

Corollary 4.8 The singularity type of any Goursat structure on a manifold of di- 
mension n belongs to the Jacquard language J„-3. Conversely, any word of J n _3 is 
the singularity type of some Goursat structure. 
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4.6 Singularity Type of the N- Trailer System 

In this Subsection, we come back to the n-trailer system, for which we compute the 



singularity type. Our study stays very close to that of Jean J29|. For the n-trailer 
system we define, following pm , the sequence of sets for i > 0, of real numbers 
by the relations 

«i = {-§ , +f } 
CKj+i = {arctansin(a), arctansin(a) + n : a G c^}. 

Note that card(«j) = 4, for i > 2. Now, consider the n-trailer system r n at a 
configuration p = (£i, £2, #0; • • • > $n) of R 2 x (S 1 ) n+1 . Define, inductively, a word 
5 T n (p) = wi ■ ■ ■ w n of J n by wi = ao and, for i = 2, . . . , n, by the relations 

u>j = 01 if 0, - 6»i_i G oti; 

iy» = a fc+1 if 6»j - 6»j_i G a k+1 and iy»_i = a fc ; (32) 
Wi = ao otherwise. 

This definition leads to a characterization of the singularity type in the coordinates 
of the n-trailer system, which coincides with the stratification of the singular locus 



given in |29[ . 



Proposition 4.9 Let V be the Goursat structure spanned by the n-trailer system 
r n on R 2 x (S 1 )^ 1 . We have 

MfO = <&r"(p)- 

Moreover, in the coordinates (£1, £2, 9q, . . . , 6 n ) of the n-trailer system, we have 

Sj ^ = {6 n _i — 6 n -i-i G Oil, . . . , n -i+j — On-i+j-l ^ a j+l}i (33) 

for < i < n — 2 and < j < i. 

Like in Section^] we will use the notation Sj = sin(0j — 0i_i) and q = cos(0j— 0j_i). 
Moreover, we define the product 7T^ by the relation 71"^ = rij'=fc c n-j+i5 if < A; < /, 
and by 7r^ = 1, if k > I. 



The proof of the Proposition [O] will use the two following Lemmas. The first 
one is analogous to Lemma |L5| , of Subsection |4.5| , it shows that the characteristic 
distributions of the n-trailer are rectified in (£i, £ 2 , 0o, . . . , n ) coordinates. Its proof 
is straightforward and left to the reader. 
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Lemma 4.10 Let T> be the distribution spanned by the n-trailer system (t^t^). 
The derived flag of T> is given by 

P(i) = (i--^.^9dr + c "-^ 1 ), forO< i <n. (34) 
The characteristic distributions Ci of T>^ +1 ^ are given by 

& = "'Kb). forO<i<n-l. (35) 

Moreover, we have 

= (&>■■■ > ^ + ^r*- 1 ). (36) 

fc=i 

/or < « < n — 1 and < j < i. 

The Lemma below is essentially a trigonometric identity and its proof, based on 
an induction argument, is also straightforward. We also leave it to the reader. 

Lemma 4.11 Let (£i, £ 2 > #o, • • ■ > #n) G IR 2 x (S' 1 ) n+1 6e a /zxed pomt o/ £/ie con- 
figuration space of the n-trailer. If 9 n - i+ k — 6 , n _ i+fc _ 1 G « fe+1 /or < k < j — 1 

i i 

^(•Sn-i+fc-i) Wi) ae n d i+k - 2 = Sn-i+j~ic n -i+j(^2 ae n _! +fc _ 2 )- 
fc=i fc=i 



Proof of Proposition |4.9| To start with, let us prove that: 



Sj — {0 n -i — &n-i-l & Oil, . . . ) @n-i+j ~ @n-i+j-l £ 

For any fixed z we will prove the result by induction on j. The relations (|34"D and 
( |55| ) imply that T>^(p) = Ci + \(p) if and only if c n _j(p) = 0. That is, if and only if 
&n-i — 6n-i-i £ which implies that the Proposition to be true for j = 0. Now, 
assume the Proposition true up to j — 1. The relation (|36| ) implies that T>^~^ is 
given by 

i+i 
fc=l 



The induction assumption, together with Lemma |4.11| , implies that for any point p 
that belongs to Sjii the linear subspace T>^ l ~^{j)) is equal to 

j 



S n -i+j d6n d x + ( S n-i+j-lCn-i+j)(^2 d9n _ d i+k _ 2 ))- 



k=l 
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The induction argument says that 



Since 



Sj-l — {@n-i — On-i-1 £ a l, • • • ? Qn-i+j-l ~ $n-i+j-2 £ Olj}. 



i+1 

T <7® — ( 9 9 9 9 9 9 9 



fe=l 

we have Z>0 _ J) (p) n TpS}!^ 7^ (p) if and only if s n _ i+i (p) = s n - i+j -i(p)c n - i+j (p). 
That is, we have p G S'j 1 '' if and only if p 6 S r ^_ 1 and 9 n ^ i+ j — Q n ~i+j-\ £ which 
ends the induction argument. 

Now, the form of obtained in the previous paragraph together with the 
definitions of 5x> and 5 r n, imply directly that 5z>(p) = S T n(p). □ 



5 Growth Vector 

5.1 Lie Flag and Growth Vector 

The Lie flag of a distribution V is the sequence of modules of vector fields V C 
Pi C ■ ■ • defined inductively by 

V = V and £> m = V { + [D , A], for i > 0. (37) 

This sequence should not be confused with the derived flag ([!]), introduced at the 
beginning of the article. In general these two sequences are different. Nevertheless, 
for any point p in the underlying manifold M, the inclusion T>i(p) C V^ l \p) holds, 
for i > 0. Note that for a Goursat structure, unlike the elements of the derived flag, 
the elements of the Lie flag are not necessarily distributions of constant rank. 

A distribution T> is completely nonholonomic if, for each point p in M, there 
exists an integer N(p) such that Pjy(p)(p) — T p M. The smallest such integer is 
called the nonholonomy degree of T> at p and we denote it by N p . For a Goursat 
structure on a manifold of dimension n, the inequality N p < 2 n ~ 3 holds for each 
point p in M (see e.g. IfffiQ. For the n-trailer system, sharper bounds were obtained 
2S], [12 1, and It follows from our Theorem |3.3j , which states that any Goursat 
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structure is locally equivalent to the n-trailer system, that they hold also for any 
Goursat structure. 
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Definition 5.1 Let P be a completely nonholonomic distribution. Put di(p) = 
dimDj(p), for < i < N p . The growth vector at p of the distribution T> is the finite 
sequence (d (p), . . . , d Np (p)). 

Recall that if at a given point a Goursat structure can be converted into Goursat 
normal form (|]) then this point is called regular and that otherwise it is called 
singular (see Section ||). The set of singular points is called the singular locus. An 
elegant characterization of this set, that emphasizes the importance of the growth 
vector in the study of Goursat structures, has been obtained by Murray flST] . A 
different characterization can be found in and [39fl. 



Theorem 5.2 (Murray) Let p be a point in a manifold M of dimension n. A 
Goursat structure on M can be converted into Goursat normal form in a small 
enough neighborhood of p if and only ifT>i(p) = T>^(p), for < i < n — 2. 



5.2 Growth Vector of the N- Trailer System 

Let d = (do, . . . , djv) be a finite sequence of integers such that do = 2, d^ = n, and 
di < d i+ i < di + 1, for < i < N — 1. The dual of the sequence d is the sequence 
d* = (djj, ... ,d*) defined by 

d* = card{j > : dj < i} + 1, for 2 < % < n. 

In other words, the integer d* indicates the first position, starting from the left, where 
the integer i appears in d. We obviously have d\ = 1 and d* = N + 1. It is trivial 
to check that each sequence d is uniquely defined by its dual d*. For example, we 
have the following dual sequences: (2, 3, 4, 5, 6)* = (1, 2, 3, 4, 5), (2, 3, 4, 5, 5, 5, 6)* = 
(1, 2, 3, 4, 7), and (2, 3, 4, 4, 5, 5, 5, 6)* = (1, 2, 3, 5, 8). 

Now, following we define a set of functions that will allow us to obtain 
a formula that gives the growth vector of an arbitrary Goursat structure at an 
arbitrary point, as a function of its singularity type at this point. We start with 
Jean's formula for the ra-trailer. Recall that J n denotes the Jacquard language 
(see Section ||) and that the shift of a word is defined by (w\ ■ ■ ■ w n )' — W\ • • • u>„_i 
and (e)' = e (we will denote [w')' by w"). 

For any i > 2, we define functions $ : U n>i _ 3 J„ — ► N . We take /^(w) = 1, 
foiw) = 2, and (3±(w) = 3, for any word w in Un>*-3*^»- If i > 5 then we define 
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inductively, for any word w in {J n>i _ 3 J n , 

AH = A-iK) + if w = (w')ai 

Pi(w) = 2 A-i(w') - A-2H') if w = HK and fc > 2; 
(3i(w) = Pi-\{w') + 1 if w = (w')a , 

For example, for the word a aiao ; we have: 

P 5 (a aia ) = f3i{a ai) + 1 = 3 + 1 = 4 

/? 6 (a aia ) = AsfaoOi) + 1 = (/3 4 (a ) + /3 3 (e)) + 1 = (3 + 2) + 1 = 6. 
An other example is given for the word for which we have: 

/3 5 (a aia 2 ) = (3 4 (a ai) + 1 = 3 + 1 = 4 

I3&{a a 1 a 2 ) = 2p 5 (a a 1 ) - /5 4 (a ) = 2(/3 4 (a ) + /3 3 (e)) -3 = 7. 
Having recalled the functions we are now able to recall the formula, obtained 



by Jean [29|, that gives the growth vector of the n-trailer system. 



Theorem 5.3 (Jean) Consider the n-trailer system at a given point p of its con- 
figuration space M? x (S 1 )^ 1 at which it has singularity type S T n(p). The sequence 
of integers (d^ip), ■ ■ ■ , d* n+z (p)) dual to the growth vector of the n-trailer system at p 
is given by d*(p) = /3i(6 T n(p)). 

5.3 Growth Vector of Goursat Structures 

The following result is fundamental. It shows that the growth vector of any Goursat 
structure is a function of its singularity type. 

Theorem 5.4 Let V be a Goursat structure on a manifold M of dimension n > 3, 
defined in a neighborhood of a given point p in M that has singularity type 8t>(jp)- 
The sequence of integers (d^ip), ■ ■ ■ , dual to its growth vector at p is given by 

d*(p)=(3 i (5 v (p)). 



Proof of Theorem 5.4 Let T> be a Goursat structure on a manifold M of dimen- 



sion n > 3, defined in a neighborhood of a given point p in M. By Theorem |3.3| . 
the Goursat structure T> is locally equivalent at p to the n-trailer system, consid- 
ered around a well chosen point q of its configuration space. By Theorem |5.3|, the 



sequence of integers (dl(q), . . . , rf* +3 (g)) dual to the growth vector of the n-trailer 
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system at q is given by d*(q) = /3i(6 T n(q)). By Proposition |Q| , the singularity type 
of the ^-trailer system at q equals 8 T n(q). Since the singularity type is invariant 
under diffeomorphisms, we have 8x>{p) = 8 T n(q). Since the growth vector is invari- 
ant under diffeomorphisms, the sequence of integers (d|(p), . . . , rf* +3 (p)) dual to the 
growth vector of V at p is given by d*(p) = Pi(8 T n(q)) = (3i(8v(p))- □ 

The latter result obviously implies the following one, which gives the formula for 
the growth vector of an arbitrary Kumpera-Ruiz normal form. 

Corollary 5.5 Let k u be a Kumpera-Ruiz normal form on R n ; for n > 3. The 

sequence (d%, . . . , d*) dual to its growth vector at zero is given by d* = /3i(8 K n). 



5.4 Growth Vector and Singularity Type 

We proved in the previous Subsection (Theorem |5.4j ) that the singularity type of 
any Goursat structure at a given point determines its growth vector at this point. 
Now, we will prove the converse of this fact. 

Theorem 5.6 Two Goursat structures have the same growth vector at a given point 
if and only if they have the same singularity type at this point. 



The proof of Theorem [OS] will be based on two Lemmas: 



Lemma 5.7 Let i and k be two integers such that i > 1 and < k < i — 1. For 
any word w in Un>i^« we have the following relations: 

(i) f3 i+i (waia 2 ■ ■ -at) = 2i + 3; 

(ii) f3 i+i {waia 2 ■ ■ ■ aj_ fe a§) = 2 i - k + 3; 
(hi) f3 i+ 4(wci ■ ■ -Ci) = i + 3; 

where Cj, for 1 < j < i, are any letters satisfying Cj ^ a\. 



Proof of Lemma |5.7| First Item. Item (i) is true if % = 1 because, for any word w 
in Un>i J n, we have 

As(wai) = /3 4 H+/3 3 «) 
= 3 + 2 = 2-1 + 3. 

It is also true if i = 2 because, for any word w in Un>i^f we nave 

/3 6 (w7aia 2 ) = 2 • /3 5 (wai) - /3 4 (w) 
= 2-5-3 = 2- 2 + 3. 
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Now proceed by induction on i > 3. Assume that Item (i) is true up to i — 1. Then, 
for any word w in {J n>1 J n , we have 

A+4<waia 2 •■•a i ) = 2 • /%_i)+ 4 (u>aia 2 • • • - /%-2)+4(waia 2 • • • a;_ 2 ) 
= 2 • (2 • (i - 1) + 3) - (2 • (i - 2) + 3) 
= 2i + 3. 

Second Item. Let us proceed by induction on i. It follows from Item (i) that, for 
i > 1, Item (ii) is true for i — 1 and k — 0. Assume that Item (ii) is true up to i — 1 
for any < k < i — 2. Then we have, for 1 < k < i — 1 and for any w in U n>1 J„, 
the following relation: 

I3 i+A (waia 2 ■ ■ ■ Oj_ fc ao) = /%_i) +4 (u>aia 2 • • • a^.i).^-!)^" 1 ) + 1 

= 2 - (i — 1) — (A; — 1) + 3 
= 2i — + 3. 

Since, by Item (i), Item (ii) is true for k — 0, it follows that Item (ii) holds for any 
% > 1 and any < fc < i — 1. 

Third Item. Item (iii) is true if % = 1 because /^(iuci) = 4 for any word w in Un>i^n 
(recall that c\ ^ ai). It is also true if « = 0. Now proceed by induction on i. Assume 
that this Item is true up to i — 1, then we have either 

A +4 (wci •••a) = 2 • /%_i )+4 (wci • • • Cj_i) - /3(j_ 2)+4 (wci • • • Cj_ 2 ) 
= 2 - ((i - 1) + 3) - ((i - 2) + 3) 
= i + 3 

or 

/3 i+4 (wci • • • q) = /3(j_i )+4 (wci • • • Cj_i) + 1 
= ((i-i) + 3) + l 
= i + 3, 

which ends the proof. □ 

Lemma 5.8 Let i be an integer such that % > 5. Consider two words W\ and w 2 of 
the Jacquard language J\, with I > % — 3 ; such that: 

(i) A(wi) > Pi(w 2 ); 

(ii) A_iK) > ft-iW); 

(hi) Pi{ Wl ) - A_iK) > AM - ft-iW). 

Tnen, /or any integer k > 1 and /or any word w snc/i i/iat u>iu> and w 2 w belong to 
J k +i, we have (3 i+k (wiw) > (3 i+k {w 2 w). 
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Proof of Lemma |5.8| Consider two words w\ and w 2 in Ji, with I > i — 3, that 
satisfy conditions (i)-(iii). Let a be any letter such that w\a and W2CL belong to Ji+x- 
Then we have the three following cases: 
If a = a then 

(3 i+1 (wia ) = (3i(wi) + l 
f3 i+1 (w 2 a ) = Pi(w 2 ) + 1. 

If a = ai then 

f3 i+1 {wiai) = pi(wi) +p i -.i(w' 1 ) 
Pi+i(w 2 ax) = Pi(w 2 ) + fti-i(w' 2 ). 

If a = cij then 

P i+ x(w 2 aj) = (3i{w 2 ) + Pi(w 2 ) - Pi-x(w 2 ). 

Therefore, in any case, the words wxa, and W2C1 satisfy the three conditions (i)- 
(iii); and in particular, we have f3 i+ i{wia) > /3 i+ x(w 2 a). Hence the Lemma is true 
for k — 1. An induction argument on the length of w, based on the same relations 
as for k — 1, ends the proof. □ 

Proof of Theorem |5.6| By Theorem |5.4| , if two Goursat structures have the same 
singularity types at p and p, respectively, then they have the same growth vector at 
p and p, respectively. Now, we will prove the converse. Suppose that w and w are 
the singularity types of two distributions T> and T> at p and p, respectively, that is 
w = Sx> (p) and w = 5x> (p) ■ We will show that if w 7^ w then there exists an integer i 
such that Pi (w) 7^ (3i (w). 

It is easy to check that if w and w are two words of the Jacquard language J n 
such that w 7^ w then there exists (after a permutation of w and w, if necessary) 
three words z, v, and v such that both w = vz and w = vz, and which satisfy either 

{v = uaia 2 ■ ■ ■ aj_fcCio 
V = UCxC 2 ---Ci, 

where < k < i — 1 and Cj 7^ ax for 1 < j < i, or 

v = ua\a 2 ■ ■ ■ aj-fcCtg 
v = ua\a 2 ■ ■ ■ di^ia l , 
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where k ^ I. 

For each of these two cases we can apply Lemma |5.7| . In the first case we have 
Pi + 4,(v) = 2i — k + 3; while l3 i+i {v) = i + 3. Since k < i—1 we have /3i + ±(v) 7^ l3 i+i {y). 
In the second case we have (3 i+4: (v) = 2 i — fc+3; while /3 i+4: (v) = 2 i — Z+3. Since k ^ I 
we have (3i + i{v) ^ /3 i+4 (5). Therefore, in both cases, we have (3i+4(v) 7^ Pi + ^(v)] but 
j3 i+3 (v) = Pi + s(v), since by the Item (iii) of Lemma [5/7| they are both equal to i + 2. 
Put i = {i + 4) + \z\. By Lemma pT%| , we have Pi (w) 7^ f3 io {w). □ 



5.5 Computing the Singularity Type 

Up to now, we have worked with a definition of the singularity type that uses the 
submanifolds Sj. Although being geometric, that is independent of a description of 
the Goursat structure in particular coordinates, it does not tell us how to compute 
this invariant (unless we know how to compute all Sj explicitly). In order to fill 
this gap we give the following Proposition, which yields to a constructive procedure 
to compute the singularity type of any Goursat structure in terms of its growth 
vector. Its proof is straightforward. 

Proposition 5.9 Let V be a Goursat structure considered in a neighborhood of a 
point p that belongs to a manifold of dimension n > 5. For < i < n — 5 and 
1 < j < h the point p belongs to Sj if and only if the growth vector at p of the 
distribution T>^~^ starts with 

(i - j + 2, i - j + 3, . . . , i + 2, i + 3, i + 4, . . . , i + 4, i + 5), 

where the integer i + 4 is repeated j + 2 times. 



6 Abnormal Curves 

6.1 Integral and Abnormal Curves 

Let M be a smooth manifold of dimension n and let A be a set-valued map A : 
M — > TM such that A(p) C T P M, for each point p in M. Note that we do not ask 
A(p) to be a linear subspace of T P M, but just a subset of T V M. Neither we ask A 
to be smooth. An integral curve of A is an absolutely continuous map x : / — > M, 
from an interval / C R to M, such that x(t) belongs to A(x(t)) for almost all t 
in /. A nontrivial lift of x(-) is an absolutely continuous map P : I — > T*M such 
that P{t) belongs to T* (t) M and P(t) ^ for each t in /. 



41 



Locally, all integral curves of a rank k distribution V — . . . , fk) can be de- 
scribed as solutions of an (underdetermined) ordinary differential equation. Indeed, 
for any given integral curve x(-) of T> we can clearly find k real- valued measurable 
functions Ui, for 1 <i < k, such that 



i(t) = J2tjMt))u i (t) (38) 

holds for almost all t in I. These functions itj are called controls. Observe that the 
controls associated to an integral curve are not uniquely defined. In control theory, 
an overdetermined differential equation of the form (|38|), where the functions Ui for 
1 < i < k can be taken as arbitrary measurable functions, is called a control system. 
Informally, the system (|3"8D can be seen as a "parametrization" of the set of all 
integral curves of T> by k real-valued measurable functions. 

Roughly speaking, a solution x(-) of fl38D is abnormal if it is a singular point 
of the end-point mapping or, equivalently, if the linearization of the control system 
along x(-) is not controllable. Many equivalent definitions of the concept of abnormal 
curves are available (see e.g. the papers [|J, 0, |65l , |f75| , the survey article ||45| , 



and the references given there). The definition that we will use is the one that 



appears in Pontryagin's Maximum principle |57|]- For further details, we refer the 
reader to the above mentioned works. 

Since the results of this section will be local we can work in a coordinate chart 
x : M — > W 1 . Denote by (x,p) the corresponding coordinates on T*M. In these 
coordinates, the Hamiltonian of the control system ( j38|) associated to a distribution 
V = (fx, . . . , f k ) is the function defined on l"xl"x R k by 

/ k 
H(x,p,u) = ( p , Y,fi( x ) u i 

\ i=l 

where both x and p belong to R™ and u = (ui, . . . ,Uk) belongs to M fc and (•, •) 
denotes the pairing between vector fields and differential forms. 

Definition 6.1 An integral curve x : I — > W 1 , corresponding to a measurable 
control u : I — > lR fc , of the control system (|38| ) is called abnormal if it admits a 
nontrivial lift (x(-),p(-)) such that 

_ dH(x(t),p(t),u(t)) 

X[) ~ dp 

= dH(x(t),p(t),u(t)) 

dx 
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and 

dH(x(t),p(t)Mt)) = 
du 

for almost all t in I. 

By definition, an integral curve of a distribution T> — (fx, . . . , /&) is abnormal if 
it is an abnormal curve of the corresponding control system. It is well known that 
the abnormal curves of T> depend neither on the choice of coordinates nor on the 
vector fields /i, . . . , fk chosen to span the distribution. 

Let / C M be an interval. For any t e / and for any e > 0, denote by I £ (to) the 
intersection I C\ [to — e,to + e\. An integral curve x : I — > M is locally abnormal if 
for each to in / there exists a small enough e > such that the restriction of x(-) to 
I e {to) is abnormal. 



6.2 Abnormal Curves of Goursat Structures 

Let P be a Goursat structure on a manifold M of dimension n > 3. Recall that its 
singularity type can be computed using the sequence of canonical manifolds defined, 
for < i < n — 5, by 

= {qeM:V^{q)=C l+l {q)} 

and, for 1 < j < i, by 

Sf = {qE Sf\ : V^\q) R T q Sf\ ± C^(q)}, 

where the distributions C« are the canonical distributions of Proposition |4.1| (see 
Section Assume that for two given non-negative integers i and j, such that 
0<i+j<n — 5 we have Sf +i) ^ 0. In this case, we can define on Sj^ a smooth 
distribution A * by taking 

A ( ;\q)=V^(q)nT q Sf +J \ 

for each point q in is easy to check, using a Kumpera-Ruiz normal form, 

that Aj is indeed a smooth distribution and that its rank is i + 1. Although each 
A^ is defined only on S^^, we can extend the definition of Aj to M by taking 
(q) = f° r an points q that do not belong to Sj and thus consider Aj as 
a set valued map defined everywhere on M. This extension allows us to define, for 
any < i < n - 5, a subset A {i) C TM by 

^ w (ff)=C(r/)u( U ^°(?) 

\ 0<j<n— i— 5 
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for each point q in M. Note that, usually, the subset A.® C TM is not a distribution. 

By definition, we take A^ n ~^ = C n -4- Moreover, we define A^ -3 ^ as the char- 
acteristic distribution of V^ n ~ 3 \ which is equal to C n _ 4 if n > 4 and equal to {0} if 
n = 3. Finally, we take A^ n ~ 2 ' = 0. Observe that the set-valued maps A^' 3 ^ 1 = {0} 
and A( n ~ 2 ^ = are different. Indeed, the first one has trivial integral curves (points); 
while the second one has no integral curves at all. 

Theorem 6.2 Consider a Goursat structure T> defined on a manifold of dimen- 
sion n and fix an integer i such that < i < n — 2. An integral curve of is 
locally abnormal if and only if it is an integral curve of A^'. 

For i = n — 4, n — 3, and n — 2, the distribution X>w is of rank n — 2, n — 1, 



and n, respectively, and the proof of Theorem |6.2| follows easily from well known 
results. Indeed, if i = n — 4 then the distribution D®, which is of rank n — 2, can 
be transformed into a direct generalization of Engel's normal form (||32[|, H3|, [73|, 



and |76|) given by Theorem |A.4j (see Appendix ^|), where we have to take k = n — 2 



and m = 2. In this case, the abnormal curves of D*™ -4 ) are clearly the integral 
curves of C n _ 4 (see Lemma |TB] below). If i = n — 3 then the distribution 
which is of rank n — 1, is annihilated locally by a 1-form w such that du A 7^ 
and (rfu;) 2 A = 0. This property is equivalent to the fact that the characteristic 
distribution of I>( n ~ 3 ) is of corank 2 in T>( n ~ 3 ) (see [|]), and it implies that p( ra_3 ) 
is locally given by the normal form of Theorem |A.4j , where k = n — 1 and m = 1. 



Note, however, that this form does not follow from Theorem [A.4| whose condition, 
when m = 1, is only necessary but not sufficient. In this case, it is straightforward 
to see that the abnormal curves of £>( n-3 ) are the integral curves of the characteristic 
distribution of T>^ n ~ 3 \ which is an involutive distribution that has rank n— 3. Finally, 
if i = n — 2 then the situation is even simpler because 

2>(n-2) = TM; which i m pii e s 
that p( n ~ 2 ) has no abnormal curves at all. Hence the only values of i that will be 



considered in the proof of Theorem are < i < n — 5 



In order to explain further the meaning of Theorem |6.2| we would like to empha- 
size the following points, relative to the geometric structure of A® and its integral 
curves. These facts follow directly from our study of the singularity type (see Sec- 
tion fD and will be used in the proof of Theorem |6.2|. 



(i) Although for each point q in M we have, by definition, 



A®(q)=Ci{q)u{ U Af(q) 

t0<j<n—i—5 
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the relations S% +k) n Sf +j) = for k ^ j (see Proposition fO| ) imply that, for a 
fixed point q, only two possibilities can occur. Indeed, we have either 



-4%) = C;(g) or A {l \q)=C t (q)UAf\q), 



for a unique integer j such that < j < n — i — 5. In other words, for each point g 
the subset A^'(q) C T g M is the union (not the sum!) of either one or two linear 
subspaces of T q M. Note that if i > 1 then Ci(q) fl A * (q) = C$_i(g). 

(ii) For < z < n — 5, we define the set Ki = {J™ZfS^; for any other value 
of i we take Ki = 0. We will call this set the singular locus of P^. If i = n — 4, 
n — 3, or n — 2 then, by definition, the singular locus is empty, which explains why 
these cases are simpler. If % = then this definition agrees with the one given in 
Section |5] for the singular locus of T>. It follows directly from the properties of the 



submanifolds Sq (see Proposition 4.4 ) that K{ is a stratified manifold. In fact, 



in Kumpera-Ruiz normal form coordinates, this set is an algebraic variety defined 
by a single polynomial equation of the form n^o^L. = 0; where the integer m 
corresponds to the number of singularities of X>W, which is, in general, smaller than 
the number of singularities of T>. For any point q that does not belong to Ki we 
clearly have A^\q) = Q(q). Note, however, that there exist points of K^ for which 
we also have A^(q) = Q(q). 

(iii) For < i < n — 5, define the set Lj = \J^Zq~ 5 Sj t+ ^ . Since for any j we have 
Sj 1 ^ C Sq \ it follows that L, C Ki. For < i < n — 5, the set of points such that 
A^(q) ^ Ci(q) is precisely Lj. In other words, the set A^'(q) is a linear subspace of 
T q M if and only if q does not belong to Lj. Unlike Ki, the set Lj is always a smooth 
submanifold of M. Note, however, that L 4 can have several connected components 
and that the dimensions of these components can be different. Nevertheless, in a 
small enough neighborhood U of any of its points, the submanifold Li is connected 
and coincides with one and only one of the smooth manifolds Sj t+ ^ fl U. 

For example, in the case of a distribution spanned by a Kumpera-Ruiz normal 
form on M. n , the set Lj is connected. If non-empty, the set is a codimension j + 1 
linear subspace of M™, where j is the only integer such that Sj t+ ^ is non-empty. In 
the case of the n-trailer system, the situation is quite different. For example, for 
the two-trailer system, the submanifold L has two connected components, given 
by {02 — &i = 7r /2} and {9 2 — 9 1 = —n/2}. Each of them has codimension 1. For 
the three-trailer system, the submanifold Lq has six connected components given, 
respectively, by {9 3 - 9 2 = vr/2}, {9 3 - 9 2 = -vr/2}, {9 3 - 9 2 = vr/4; 9 2 -9 x = vr/2}, 
{9 3 -9 2 = -3tt/4; 9 2 -9 x = vr/2}, {9 3 -9 2 = -vr/4; 9 2 -9 x = -vr/2}, {9 3 -9 2 = 3vr/4; 
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9 2 —9i = 



— 7r/2}. Two of them have codimension 1; four of them have codimension 2. 



We consider now a more detailed example. Let T> be the distribution spanned 
by the following Kumpera-Ruiz normal form on R 7 : 

where c-j is either equal to or 1. When c 7 = 1 the singularity type of T> at zero is 
and the growth vector at zero is (2, 3, 4, 5, 5, 6, 6, 6, 7); while when c 7 = 
the singularity type is aoOiOifl^ and the growth vector (2, 3, 4, 5, 5, 5, 6, 6, 6, 6, 7). 
In both cases, we have 

= 0, Sf = {x 6 = 0}, and = {x 5 = 0}. 

Therefore, the singular loci of the distributions T>^\ PW, and are given re- 
spectively by 

K = {x 6 x 5 = 0}, K x = {x 6 x 5 = 0}, and K 2 = {x 5 = 0}. 

If C7 = 1 then, in a small enough neighborhood of zero, we have S[ 1] = 0; but if 
c 7 = then we have Si = {x-j = xq = 0}. In both cases we have = 0. 

If C7 = 1 then, in a small enough neighborhood U of zero, we have A^ = Cq = 
(9^ )' wmcn is a smooth distribution on U; but if c 7 = then the subset A^ 
coincides with the smooth distribution C = (gf-) outside L = {x 7 = x e = 0} while 
for any point p of Lq we have 

^ (0) (p) = (^)(p)u(^)(p), 

which is clearly not a distribution. In both cases, we have A^ = C\ = (7^, 7^) 
outside Li = {xq = 0} while for any point p of L\ we have 



Finally, we have A^ = (^-, gf^, ^) outside L 2 = {x 5 = 0} while for any point p 
A \p) = {g^, q^, d^)(p) U a^)(p)- 



> 8x7 ' (9X6 ' &5 

of L 2 we have 



We proceed now to the proof of Theorem |6.2| , which states that an integral curve 
of PW i s locally abnormal if and only if it is an integral curve of A^\ The proof 
will be based on the three following Lemmas. 
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Lemma 6.3 Consider a Goursat structure T> defined on a manifold of dimension n, 
and fix an integer i such that < i < n — 4. An integral curve o/PW that has an 
empty intersection with the singular locus Ki is locally abnormal if and only if it is 
an integral curve of Ci, and thus of . 



Proof of Lemma B.3| Let 7 : I — > M be an integral curve of that does 



not intersect the singular locus K^. Since we are outside Ki it is easy to show, 
using a direct generalization of Goursat 's normal form (see and P5|), given 
by Theorem |A.4j, that for any fixed to in I we can find a local coordinate chart 



x : U — > R n centered at 7 (to) an d such that: 

T)W — I d 9 ™ 9 1 1 ™ 9 1 9 

^ 5o:i ' Xl+1 ax l+2 1 *~ Xn ~ 2 dx n -! ^~ dx n 

Chose a small enough e > such that the restriction of 7 to I e (to) is completely 
contained in the open set U. Then, the curve x o 7 : I e (t ) — > R™, which will be 
denoted shortly by x(-), is almost everywhere a solution of the following control 

system 

X\ = Ui 



X{ Ui 

x i+1 = u i+1 m 

%i+2 — %i+l u i+2 



Xn-l — X n -2Ui + 2 

Since the coordinate chart is centered at 7(^0); we have x(to) = 0. The Hamiltonian 
of this system is given by 

H(x,p,u) = Y!k=\PkU k + J2k=i+2Pk X k~lUi+2 + PnUi+2- 

Therefore, the curve x(-) is abnormal if and only if there exists a non-trivial lift 
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(x(-),p(-)) that satisfies, almost everywhere, the following differential equation 



Pi 







Pi 







Pi+l 



-p i+2 U i+2 



(40) 



Pn-2 = -Pn-lUi+2 
Pn-l = 
Pn =0 



and, moreover, is such that pk — 0, for 1 < k < i + 1, and p n = — J2t=i+2Pk x k-i- 
The latter condition is a consequence of |^ = 0. 

Necessity. Assume that x(-) is not an integral curve of C«. We will prove that x(-) 



is not an integral curve of C iy there exists a measurable subset I C / E (to) such that 
the Lebesgue measure of Iq is not zero and u i+2 {t) ^ for each t in I . If x(-) is 
abnormal then p(-) is such that pj+i(t) = for each t in I £ (to). Therefore, we have 
= almost everywhere on I . Indeed, note that if an absolutely continuous 
function / on I is such that f(t) = for almost all t in I then f'(t) = for almost 
all t in J . But pi + \ = —p i+2 Ui + 2 and u i+2 ^ imply p i+2 = almost everywhere 
on Jo, which gives pi +2 = almost everywhere on Iq. We can repeat the previous 
argument to obtain p^ = 0, for 1 < k < n — 1, almost everywhere on J . Since 
p n = — J2k=l+2Pk x k-i> we have also p n = almost everywhere on I . This gives 
Pk = 0, almost everywhere on 7 ! f° r 1 < k < n, which is impossible since p must 
be non-trivial. 

Sufficiency. Now, assume that x(-) is an integral curve of Cj. In order to prove 
that x(-) is abnormal, we will consider the lift defined by pt = for 1 < k < n, 
with the exception of for which any non-zero real constant can be taken. 

Since x(-) is an integral curve of Cj we must have Ui +2 (t) = almost everywhere 
on I e (to), which implies that p(-) satisfies (0). Moreover, since x(to) = 0, we have 
Xk{t) = Xk{to) = 0, for each t in I £ (to) and for i + 2 < k < n. Thus p n , which 
was taken to be zero, satisfies p n = —^2k=i+2Pk x k-i (recall that pi +2 = 0). In other 
words p(-) satisfies ^ = 0. Finally, since p n -\ ^ 0, our lift is non-trivial, which 
implies that x(-) is abnormal. □ 



is not abnormal. In the coordinates of (^) we have = ( 



a a 



). Since x(-) 



dxi ' ' ' ' ' dxi+i 
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Lemma 6.4 Consider a Goursat structure T> defined on a manifold of dimension n 
and fix an integer i such that < i < n — 5. Let x(-) be the restriction of an integral 
curve ofT>^ to the interval I e (to); where e > 0. If a fixed measurable subset / C 1 
is such that for each t in I H / e (to) we have x(t) A^\x(t)) then, for a small 
enough e > 0, we have x(t) ^ Ki for almost all t in IqD I £ (to). 



Proof of Lemma |6.4j Let x : I £ (to) — > W 1 be the restriction to the interval I £ (to) 



where e > 0, of an integral curve of D". Firstly, if x(t ) ^ Ki then there exists 
a small enough e such that x(-) does not intersect Ki and thus, in this case, the 
Lemma is trivially true. Secondly, if the Lebesgue measure of Iq is then the Lemma 
is also trivially true. Finally, if the closure of Iq does not contain to then for a small 
enough e the Lebesgue measure of Iq fl I £ {to) will be and thus the Lemma will be, 
once more, trivially true. Hence, from now on, we will only consider curves such 
that x(t ) belongs to K iy the Lebesgue measure of I is not 0, and the closure of I 
contains to- Moreover, once a small enough e > has been fixed, we will denote 
also by I the intersection I n I e (to)- That is, we will assume that 7 C I £ (t ). 

For any such integral curve x(-) of it is easy to prove, using a direct gen- 
eralization of Kumpera- Ruiz's normal form, given by Theorem [A.5| (with a double 



indexation of coordinates, like in Corollary |2.4| ), that there exist coordinates on IR n 
in which x(-) is a solution of the following control system: 



x\ = U\ 

x\ = u 2 



#1+1 — ^2 + 1 



(41) 



for i + 2<g</c ifp = and 
o<r<p-\ Kr ) ' q q ' " for 2 < q < k p if 1 < p < m 



±P q=\ EI 4 r ( X g-1 + 



=| EI x k r I u *+2 fo r < p < m, 

0<r<p-l 

where x = (x^x^,... , x^ ,x\,... , x™ , x™ +1 ) and x(t ) = (recall that we 
assume that x(t ) belongs to the singular locus, which means that m > 1). Moreover, 
the integers ki that appear in (|4T| ) satisfy i + l<A; <n — 4 and fci > 1, . . . , k m ^i > 
l,k m > 3,k m+ i = 1 and YlT=o^r = n. Observe that the number m > 1 is the 
number of singularities of T>^\ which can be smaller than the number of singularities 
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of V. In these coordinates, the singular locus is given by 

If c° q = 0, for all % + 1 < q < ko — 1, then the only integer j such that <Sj contains 
zero is 7 = ko — (i + 1). Thus Li is given (see Proposition |4.4|) by 

r. _ cr(*D-l) _ f_0 _ r _...__0_ n l 
^ — D fc -l-i ~~ \ X i+l ~ X i+2 — — x k ~ U J • 

Note that if for some i + 1 < q < ko — 1 we have c° q 7^ then the submanifold Lj 
does not contain zero (it is locally empty at zero). 
For each point p of M n , we have 

C(P) = (^,-,5^)(P)- 
If c° = 0, for z + 1 < q < ko — 1, then we have 

^ {1) (P) = • • • > 5^)(P) U • • • , £r)(p), 

for each point p in Lj and ,4y'(p) = Cj(p), outside Lj. If for some i + l < q < k — 1 we 
have c° then, in a small enough neighborhood [/ of zero, we have.A (i) (p) =Ci{p) 
for each point p in U. 

Recall that the integral curve x : I e {to) —> ^ n of PW is such that x(t ) = 0. 
Assume that, at a given t of I e {to), the velocity exists and satisfies (^TJ). Then, 
if e is small enough, the velocity x(t) belongs to A^(x(t)) if and only if we have 
M i+2 (t) = or the three following conditions hold: (i) m+i(t) = and (ii) x q (t) = 0, 
for i + 1 < q < k , and (iii) c° = 0, for i + 1 < q < ko — 1. 

Now, suppose that for each £ in 7 C I £ (to) the velocity i;(t) exists and is such 
that x(t) <£ A (i) (x(t)). Recall that we can assume that the Lebessrue measure of 
Io is not and the closure of io contains to- For each t in Jo we have -u.j +2 (t) 7^ 0. 
Additionally: (a) If c° = for i + l<q<k — 1 then we can represent the subset 
Io as Io — I\ U h (with, in general, a non empty intersection of I\ and J 2 ), where Ii 
is the set of points where (i) is not satisfied and 7 2 is the set of points where (ii) is 
not satisfied, (b) If there exist an integer i + l<g<fco — 1 such that c° q 7^ then 
Io = {t G I £ {to) : Mj+ 2 (t) 7^ 0}, provided that e is small enough. We are going to 
show that, in both cases, we have x(t) ^ Ki, for almost all t in Jo- 
Case (a): Subset I±. For each t in I\ we have both w i+2 (t) 7^ and Wj+i(t) 7^ 0. 
Therefore, we have x° i+l 7^ almost everywhere on Ji, which implies that 7^ c° +1 
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almost everywhere on I\. Indeed, note that if an absolutely continuous function / 

on 1\ is such that fit) 7^ for almost all t in I\ then, for any constant c, the measure 

of the set {t E I\ : fit) = c} is zero. 

Now, using an induction argument we will show, successively, that 7^ c° +1 , 
™o _/ „o ™0 _/ „o ™0 ^ n „i ^ „i ™i J- A t 1 =L n 

£^7* 7^ 0, almost everywhere on I\. Suppose that this assumption is true up to 
Xq-i- We have two cases: either q < k p or q — k p + 1. If g < fc p then x v q = 
n <r<p-i(4 r )(^-i + Since x r kr ^ 0, for < r < p - 1, and a^_ x ^ 

and Mi + 2 7^ 0, almost everywhere on Ji, we have x p q {t) 7^ for almost all t in I\. This 
implies, almost everywhere on J 1; that 7^ if g < fc p — 1 or that x p q 7^ if q = k p . 
Otherwise q = k p + 1 and in this case x\ +l = Ylo<r<p-i( x kr) u i+ 2 - Since x r kr 7^ 0, 
for < r < p — 1, and u i+2 7^ 0, almost everywhere on I x , we have 7^ for 

almost all t in Ji, which implies 7^ c^ +1 almost everywhere on I\. This ends the 
induction argument. In particular, we have proved that x r kr (t) 7^ for almost all t 
in Ii, for each < r < m — 1. Now, recall that the singular locus is given by the 
relation n^o lx fc r = 0- ^ thus follows that we have x{t) (fc Ki for almost all t in I\. 

Case (a): Subset I 2 . We can represent the subset I 2 as I 2 = I 2 +1 U • ■ ■ U I 2 , where 
I 2 = {t E Io '■ Xq(t) 7^ 0}. Observe that, in general, the intersection of these subsets 
will be non empty. Now on each subset J|, of positive Lebesgue measure, we can 
follow the same proof as for the subset I\, starting the induction argument with x q . 
For each one of these subsets the conclusion is the same: we have x(t) ^ Ki for 
almost all t ui I 2 . 

Case (b): We have u i+2 (t) 7^ for each t in I and, moreover, there exists q such 
that c q 7^ 0, where i + 1 < q < ko — 1. Since c q 7^ 0, we can take a smaller e > 0, 
if necessary, in order to have x q+1 (t) = (x q (t) + c q )ui +2 7^ for each t in J . This 
implies that x q+l 7^ c q+l almost everywhere on J . The rest of the proof follows 
like for the subset I\ and the conclusion is the same: we have x{t) ^ Ki for almost 
all t in Io- The only difference being that the induction argument starts with x q+1 
instead of x® +l . □ 

Lemma 6.5 Consider a Goursat structure T> defined on a manifold of dimension n 
and fix an integer i such that < i < n — 5. An integral curve of that intersects 
the singular locus Ki is locally, in a small enough neighborhood of any point of 
intersection, an abnormal curve ofT>^\ 



Proof of Lemma |6.5| Let x : I £ (to) — > M. n be the restriction to the interval I £ (to), 



where e > 0, of an integral curve of that intersects the singular locus at x(t ). 
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It is easy to prove (see Lemma |B.2| ) that, for a small enough e > 0, there exist 
coordinates such that the integral curve x(-) is a solution of the following control 
system: 

X\ = U\ 

Xi Ui 

x i+l = u i+l 

x i+2 = + Ci + \)Ui + 2 

'■■ (42) 
x k = (^fco-l + c fco _i)u i+ 2 

x k +l = x k Ui+2 

x k +2 = Ui+2 

x k +3 = x k x k Q +2 u i+2 

x k +4 = x k x k Q +3 u i+2 

±j = x ko (f)j(x ko+3 )u i+ 2 for k + 5 < j < n, 

where x = {xi,x 2 , . . . , x n ) and x fco+3 = (x fco+3 , x fco+4 , ... ,x n ). Since x(t ) belongs 
to Ki, we can assume that x(t ) = 0. Moreover, like in the proof of the previous 
Lemma, we have i + l<ko<n — 4. 

The Hamiltonian of this system is given by 

H(x,p,u) = jy^iPjUj + Yl k j°=i+2Pj( X j-l + c j-l) u i+2 + Pk +l x k Ui+2 + Pk +2U i+ 2 
+ (Pk +3 x k +2+Pk +4 x k +3 + E"=fc +5PA(^o+3)J x k U i+ 2, 

which implies that any abnormal lift (x(-),p(-)) of x(-) must satisfy 

pj = for 1 < j < i (43) 

Pj = -Pj+i u i+2 for i + 1 < j < k - 1 

Pk = y-Pko+1 - Pk +3 x k +2 ~ Pko+iXko+3 ~ E"=fc +5Pi^(^o+3)) U l+2 
Pko+1 = 

Pj = -^(^fe,)+3,p)^fco M i+2 for h + 2 < j < n, 

where the *0/s are some functions of x ko+3 , x feo+4 , ... ,x n and Pk +3,Pk +4, ■■■ ,Pn, 
for ko + 2 < j < n. The exact form of these functions is irrelevant for our purpose. 
Any abnormal lift (x(-),p(-)) of x(-) must also satisfy the relation ^ = 0, which 



52 



implies pj = 0, for 1 < j < i + 1, and 

Pk +2 = -Y.%i+2PA X i-± + C i~l) + 

- (Pfco+l + Pfe +3^fco+2 + Pko+AXko+3 + E"=fc +5.Pi ( M^o+3)) ^fc - 

Recall that x(-) is an integral curve of X>W. Like in the proof of Lemma |b.4] , if 
for a given i the velocity x(t) exists and satisfies fl42|) then we have G .4^(:r(i)) 
if and only if u i+2 (i) = or the three following conditions hold: (i) u i+ i(t) = and 
(ii) Xj(t) = 0, for i + 1 < j < fc , and (iii) = 0, for i + 1 < j < k — 1. Now, 
assume that x(-) is an integral curve of .4^, that is i;(t) G ^4^(x(t)) for almost all t 
in I £ (to). In order to prove that x(-) is abnormal we must construct a non-trivial 
abnormal lift (x(-),p(-)) of x(-). Take = 0, for 1 < j ' < n, with the exception 
of Pk +A, for which we take any non-zero real constant. It is straightforward to check 
that our lift satisfies (|43|) . Indeed, the coordinate Xfc 0+3 is constant because we have 
Xfc 0+3 = Xk Xk +2Ui+2', and u i+2 (t) = or Xk (t) = for almost all t. Moreover, since 
x(0) = 0, we have Xk +3(t) = 0, for each t in I £ (to). It is also trivial to check that 
our lift satisfies ^ = 0. Since Pk +A ^ our lift in non-trivial. It follows that the 
integral curve x(-) is abnormal. □ 



Proof of Theorem |6.2| Let x : I £ (to) ~^ ^ n be the restriction to the interval 
J e (to) of an integral curve of For e > small enough, we can apply both 

Lemma |6.3| and Lemma |6.5| , which imply that if the curve x(-) is such that x(t) 
belongs to A^(x(t)) for almost all t in I e (to) then x(-) is abnormal. In other words, 
the integral curves of are locally abnormal. 

Now assume that, for a fixed interval I e (to), the curve x(-) is abnormal. Define 
the subset I C I e (t ) by 

I = {te I £ (t ) : x{t) exists and x(t) £ A {i) (x(t)} 

We will show that if e is small enough then the Lebesgue measure of Iq is zero. We 
can decompose Jo into Jo = I\ U I2, where 

Ji = {t G J : x(t) G Ki] and J 2 = {t G J : x(t) ^ J^}. 

On the one hand, the measure of I\ is equal to zero. Indeed, since for each t in J we 
have x(t) ^ 4. ( -^(x(t)), for e small enough Lemma implies that we have x(t) ^ J\"j 
for almost all t in Jo. But, on the other hand, the measure of J 2 is also equal to 
zero. To see this, let us write J 2 as J 2 = Jo fl J3, where 

J 3 = {£ G J 6 (t ) : ^ JCJ 
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Since Ki is closed (see the discussion following the statement of Theorem |6.2| ) and 
x(-) is continuous, it is clear that we can decompose J3 into a union of disjoint 
open intervals I3 = [JJ a such that, on each of them, the curve x(-) has an empty 
intersection with the singular locus K^. Moreover, since the set ^3 is an open subset 
of R the union can be taken to be countable. Now Lemma |6.3| implies that, for 
each a, we have x(t) E A^ l \x{t)) for almost all t in J a because x(-) is abnormal 
and we are outside the singular locus. Hence, since the measure of I2 is the sum of 
the measures of the sets Iq D J a (the union is countable) and the measure of each of 
these sets is zero, the measure of I2 equals zero. □ 



6.3 Abnormal Curves and Singularity Type 

Theorem 6.6 Let T> and T> be two Goursat structures defined respectively on two 
manifolds M and M , both of dimension n. Fix two points p and p of M and M, 
respectively. There exists a diffeomorphism if, with p = ip(p), between two small 
enough neighborhoods of p and p that transforms, for < i < n — 4, the abnormal 
curves of T>^ into the abnormal curves of T>^ if and only if the singularity type 
of T> at p equals the singularity type ofT> at p. 



Proof of Theorem |6.6| Necessity: Consider two distributions T> and V, defined on 



two manifolds M and M, respectively, that have different singularity types w and w 
at p and p, respectively, that is w = St>(p) and w = S-gip). We have already pointed 
out (see the Proof of Theorem |5.6| ) that if w and w are two words of the Jacquard 
language J n such that w 7^ w then there exists (after a permutation of w and w, if 
necessary) three words z, v, and v such that both w = v z and w = vz, and which 
satisfy either 

v = ua\a2 ■ ■ ■ ai-fc^o 

v = UCiC 2 '--Ci, 

where < k < i — 1 and Cj 7^ a± for 1 < j < i, or 

v = ua\a2 ■ ■ ■ a^aQ 
v = tta 1 a 2 ' ' ' di-i^Q, 

where k < I. 

In both cases, consider the abnormal curves of 2)(* 0+fc ) and T>^ Q+k \ where i$ = \z\. 
It follows directly from the definition of the singularity type (see Definition fO|) 
that for V we have p G S^°^\ ^ while for V the point p does not belong to any 
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submanifold £• 0+ . Therefore, the subset A^ 0+k \p) is not a linear subspace of 
T P M while the subset A {io+k) {p) is a linear subspace of TpM. For each vector r p of 
A^ 0+k \p) there exist an abnormal curve of p(* 0+fc ) that is tangent to r p ; for each 
vector f p of *4.(* 0+A: ) (p) there exist an abnormal curve of 2)(* 0+fc ) that is tangent to f p . 
It follows that no diffeomorphism can transform the abnormal curves 

of p(*>+*) i n t 

the abnormal curves of T>^ 0+k \ locally at p and p. 



Sufficiency: Now, assume that the singularity type 5x>(p) of V at p and <%(j5) 
of X> at p coincide and are equal to w. The distribution V (respectively T>) can 
be converted into a Kumpera-Ruiz normal form K n (respectively k n ) centered at p 
(respectively p) via a diffeomorphisms (respectively 0). Let 
(respectively x = (xi, . . . , x n )) denote the coordinates in which K n (respectively k 
is expressed. By Corollary [4.6| and the invariance of the singularity type we have 
8n n — Sa» = w. Moreover, by Proposition |4.4| , a submanifold Sj contains zero 
if and only if the submanifold Sj contains zero, which is the case if and only 
if w — W\ai ■ ■ ■ a j+1 W2, for some words W\ and w 2 such that |to 2 | = i. If those 



j J 
n \ 



manifolds contain zero then, once again by Proposition [4.4|, they are respectively 
given by 

Sj ^ = {x n —i—j = 0, . . . , x n —i = 0} and Sj' ^ = {x n —i—j = 0, . . . , x n —i = 0}. 

Now, for each integer i, we must distinguish two cases. First case: If for each in- 
teger j the submanifolds Sj and S^ + ^ are empty, in a small enough neighborhood 
of zero, then, by Theorem |6.2| , the abnormal curves of (respectively P") are, 
in a small enough neighborhood of zero, the integral curves of d (respectively Ci). 
Moreover, we have 

d d \ „„j n — / d 



d = ^- , . . . , ^- and C i: 



Oxn ' ' dx n —i J y 0x n ' ' dx n —i 

Second case: If for some integer j the submanifolds S^^ and S'j i+J, ' ) contain zero 



then, by Proposition \i.'2\ , this integer j is unique. By Theorem the abnormal 
curves of Z>W (respectively X>M) are, in a small enough neighborhood of zero, the 
integral curves of A^ 1 (respectively Jffl)- Moreover, we have 

=(&;•■■■• j£- ) (?) ^'(5) = ( jfc . ■ ■ ■ . sb ) «>. 

for each point g (respectively g) that does not belong to Sj (respectively Sj l+ ^), 
and 

■4j (?) = ( af^ > • • • > fe~7T ' fen-'!-,-! ) (?) u ( af^ ' • • • ' ) (?) 
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•Aj\o) - ( dl n ' ■ ■ ■ ' aeJL+i ' &z„-!_j-i ) (^) u ( &L ' • • • ' ) (^)) 

for each point q (respectively q) that belongs to S^ + ^ (respectively S^ + ^). 

Let $ be the local diffeomorphism of IR n defined by Xi = Xi, for 1 < i < n. In 
both cases, the diffeomorphism 

(f = (j)^ 1 O $ O 

transforms the integral curves of A^ into the integral curves of AS 1 ', and thus, by 
Theorem |6.2| , the abnormal curves of T>^ into the abnormal curves of pW. □ 



6.4 Rigid Curves of Goursat Structures 

The concept of rigidity for integral curves of distributions was introduced by Bryant 
and Hsu |7| . Rigid curves are always abnormal but there exist abnormal curves that 
are not rigid (see e.g. 0,0, and [[75]). Nevertheless, we will prove that in the case 
of Goursat structures these two concepts coincide (for C l immersed curves). 

Definition 6.7 Let Dbea completely nonholonomic distribution defined on a man- 
ifold M. Fix a closed interval [a, b] and two points p and q in M. Denote by O p>q the 
space of all C 1 integral curves x : [a, b] — > M of T> such that x(a) = p and x(b) = q, 
endowed with the C 1 -topology. An integral curve x(-) that belongs to O p>q is rigid 
if there exists a small enough neighborhood V of x(-) in O p , q such that any curve 
x : [a, b] — > M contained in V is a reparametrization of x(-). 

Roughly speaking, a curve x : [a, b] — > M is rigid if it is an isolated point of 
@x(a),x{b)- Our study of abnormal curves leads easily to the following result, which 
characterizes immersed rigid curves. This result gives also, for Goursat structures, 
a more intuitive view of the concept of abnormal curve. 

Let / C 1 be a closed interval. For any t G / and for any e > 0, denote by 
I e (t ) the intersection I H [to — £> t + e]. An integral curve x : / — » M of £> is locally 
rigid if for each to hi ^ there exists a small enough e > such that the restriction 
of x(-) to I £ (to) is rigid. 

Theorem 6.8 Let x(-) be a C 1 immersed integral curve of a Goursat structure T>, 
defined on a manifold of dimension n. The three following conditions are equivalent: 

(i) The curve x(-) is locally abnormal; 

(ii) The curve xi-) is locally rigid; 
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(iii) The curve x(-) is either an integral curve of Cq or an integral curve of A^_ lt 
for some 1 < k$ < n — 4. 

We supposed in this Theorem that the integral curve is immersed, which means 
that its velocity (defined everywhere, since the curve is C 1 ) never vanish. This 
assumption is fundamental. Indeed, an immersed rigid curve can loose its rigidity if 
we change its parametrization in such a way that it is not immersed anymore (see 
e.g. jnj). Observe also that the Theorem is stated for integral curves of T>^ but 
not for those of T>^\ if % > 1. In fact, the abnormal curves of such that their 
velocity does not belong to Cj_i have only a weaker form of rigidity: all curves that 
are close enough to them in the C 1 topology stay in a submanifold of the original 
manifold. We will consider this situation in a forthcoming work. 



Our proof of Theorem |6.8| is mainly based on the ideas introduced by Bryant 



and Hsu |7]] and Zhitomirskh [[75] . In particular, it is a direct consequence of Zhito- 



mirskii's work that the immersed integral curves of Co are rigid. 

To prove the rigidity of the integral curves of *4.&o_i we f° uow the main ideas 



of [75]. Note, however, that the statement for -4.[°' ) _ 1 is not implied by any of the 



results of [0, 0, j65|, or [JfjJ because Goursat structures are highly non-generic and 
do not fit into the large categories of (generic) rank two distributions studied in 
those papers. We would like to point out that, in the particular case of dimension 
five, the rigidity of the immersed integral curves was already observed in f48| 



Moreover, the equivalence of Items (ii) and (iiii) of Theorem |6.8| has already been 
announced in |54|| . 

Our proof of Theorem |6.8| will use the following Lemma, which will be proved 
later in Appendix [B[ The normal form that we introduce in it is analogous to the 
one used in []75J to prove that the integral curves of Co are rigid. 



Lemma 6.9 Let D be a Goursat structure on a manifold M of dimension n > 5. If 
the singularity type ofD at p is equal to 5x>{p) = wax ■ ■ ■ ak for some 1 < k < n — 4, 
where w is an arbitrary word of J n -k -i, then T> is locally equivalent at p to the 
distribution spanned by a pair of vector fields that has the following form: 

P - JL 

?1 d yi 

n 

6 = Viq^ + ■■■ + Vkog^ + + ^Vko+id^; + E <Pi(v)w 

where the coordinates yi, ■ ■ ■ ,y n are centered at p. In these coordinates, the canonical 
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submanifold Sj*°_i is given by 



S { k k o °_- 1 1) = {y 1 =0,...,y ko = 0}. 

Moreover, we have Co = (£i), for any point p ofM. n , and Aj^_i(p) = (^2)(p) ; for any 
? (fco 



pointp ofSF?*'. 



Proof of Theorem |6.8| It is well known that rigidity implies abnormality (see [0], 
|7[], and ]75| ) and thus, that (ii) implies (i). By Theorem |6.2| , any abnormal curve 
of £><°) is an integral curve of ^t (0) . Recall that A {0) {p) = C (p) U ^_i(p), for a 
unique 1 < k < n — 4, and that Co(p) fl ^°_i(p) = 0. Therefore any C l immersed 
abnormal curve of T>^ is either an integral curve of Cq or an integral curve of *4-j^_i- 
Hence (i) implies (iii). 

What remains to prove is that if a C l immersed integral curve y : I e (io) ~~ > ^ n of 

£>(0) 

is an integral curve of either Cq or A^_ 1 then it is rigid. This result is known |75|] 



for the integral curves of Cq. We can thus assume that y(-) is an immersed integral 
curve of A^_ ± (which then, by definition, stays in S^ o °_i ). It follows from Lemma 
5T9| that we can find coordinates such that y(-) satisfies yit^ — e) = and is a solution 



of the following control system: 

yi = ui 
m = y\u 2 

2/fco+i = yk u 2 (44) 

m +3 = \yl 0+ iU2 

y% = Vi(y) u 2 for k + 4 < i < n, 

with Ui(t) = (because yi(t) = on S^°~^ ) and u 2 (t) ^ for each t in I £ (t ) 
(because the curve is immersed). Since the coordinates of Lemma |T9] are chosen to 
be centered at y(t — e) = 0, from yk {t) = we conclude that y& +i(£) = 0, and 
thus that £/fc 0+3 (t) = 0, for each t in I e (to). 

Now, consider a C 1 immersed integral curve y : I e (t ) -> M™ of X>(°) that has 
the same end-points as the curve y(-). In particular, we have yk +3(to — &) = and 
yk a +3(to + e) = 0. By taking a small enough neighborhood of y(-) in O y{ t ( ^ £ ) 4ji t 0+£) 
(which is not the same as taking a smaller e > 0), we can assume that u 2 {t) 7^ 
for each t in I £ (t ), where U\ and u 2 denote the controls for which y(-) is a solution 
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of Without loss of generality, we can assume that u^it) > (the proof for 

v,2(t) < is identical). Since we have \y% +\Uz{t) > for each t in I £ (t ) and both 
Vk +3{to - e) = and yfe 0+3 (*o + e) = 0, we must have yl 0+1 {t) = 0, for each t in 
I £ (to). Together with «2 > 0, the latter relation implies that, for 1 < i < ko + 1, we 
have = for each t in I £ {to), which clearly implies U\{t) = for each t in I e (to). 
Hence, the curve y(-) is a reparametrization of the original curve y(-). Indeed, these 
two curves are C 1 immersed integral curves of £2 and have the same end points (see 



75] for more details about this last point). □ 



6.5 Rigid Curves of the N-Trailer System 

Let us illustrate Theorem by applying it to the n-trailer system. Let T> be the 
Goursat structure spanned by the n-trailer system r n on IR 2 x [S 1 ) n+1 . Recall that 
Jean's sequence of sets of real numbers aij, for i > 0, is defined by the relations 

«i = {-!>+!) 

a i+ i = {arctansin(a), arctan sin (a) + 7r : a G Q;}. 



By Proposition [0], we have 



Sj ^ — {0 n -i — 6 n -i-l G Q!l, . . ■ , On-i+j — On-i+j-l ^ 

for < i < n — 2 and < j < i (recall that n is the number of trailers, not the 
dimension of the configuration space!). It obviously follows that 

Sj = {0 n -j — On-j-l £ «1, • • ■ j 9n — On-1 £ a j+l}i 

for < j < n — 2. Each submanifold Sj has clearly codimension j + 1. Moreover, 
these manifolds are pairwise disjoint. Thus a given point p is either in none of the 
submanifold Sj at all or in one and only one of them. 

Recall also that Af\p) = P«(p) n T p Sf and that A®(p) = C {p) U Af\p), 
for a unique < j < n — 2. The canonical line field Co is given on IR 2 x (5' 1 ) n+1 
by (gf-). A simple computation shows that, on each submanifold Sj, the line field 
Af i^given by Af = (^ + • • • + ^). 

By Theorem |6.8| , a C 1 motion of the n-trailer for which the velocity never vanishes 
is rigid if and only if (i) it is an integral curve of Cq or (ii) it is an integral curve of 
A^. In the second case, the motion lies in Sj . In fact, there is an easy way to 
visualize these rigid trajectories: 
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Corollary 6.10 An immersed motion of the n-trailer system is locally rigid if and 
only if it fixes the positions in the (£1, £2) -plane of the centers of the axles of at least 
two trailers. 

For example, there passes through any configuration of IR 2 x (S l ) n+1 an integral 
curve of Cq. The corresponding motion fixes the positions in the (£1, ^-plane of all 
trailers (we just turn the front wheels). If a configuration is such that 9 n — n _i e ol\ 
(it belongs to Sq ) then, besides the motions associated to C , there is an additional 
motion given by for which the positions in the (£1, ^-plane of all trailers, 
excepted the first one, are fixed. For these motions, the center of the first trailer 
moves on a circle around the center of the second trailer, which turns with its center 
fixed (see e.g. Figure §). Observe that such a motion is possible if and only if 

9 n - n _! = ±7T/2. 



7 Contact Transformations 

7.1 A Singular Version of Backlund's Theorem 

Let T> and T> be two Goursat structures defined on two manifolds M and M, respec- 
tively, of dimension n > 3. A (generalized) contact transformation (of order n — 2) 
is a smooth diffeomorphism between M and M such that (</>*P)(p) = T>(p), for 
each point p in M. Such transformations are called automorphisms in the work of 
Kumpera and Ruiz |32| (see also [19 and [0). In a neighborhood of a regular point 



our definition coincides with the classical definition of a contact transformation on 
the space J n_2 (R, R) of (n — 2)-jets of functions that have one dependent and one 
independent variable (see || and From now on, unless we want to distinguish 

generalized contact transformations from the classical ones, we will omit the word 
"generalized" . 

Fix two points p and p of M and M, respectively. Let be a local contact 
transformation between T> and D such that <fi(p) = p. Fix two small enough 
neighborhoods U and U of p and p, respectively, such that U = <fi{U) and such 
that T> on U and T> on U are equivalent to two Kumpera-Ruiz normal forms K n 
and k n centered at p and p, respectively, and defined on two open subsets x(U) 
and x(U) of M. n , where x and x denote coordinates that transform the Goursat 
structures T> and T> into their Kumpera-Ruiz normal forms K n and R n , respectively. 
We can assume, without loss of generality, that the first prolongation (in the se- 
quence of prolongations that define n n and k n ) is regular. Namely k 4 = Rq(k 3 ) 
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and k A = Ro(k 3 ). Once such a pair of Kumpera-Ruiz charts (x,U) and (x,U) 
has been fixed, we can associate to the contact transformation a unique contact 
transformation $, between k u on x(U) and k n on x(U), by taking 

$ = 5 o o a; -1 . 

In other words x = ($oi)o^, where -0 denotes the inverse of the diffeomorphism 0. 
Observe that, since the Kumpera-Ruiz charts x and x are centered at p and p, 
respectively, we have $(0) = 0. We will denote by $j the z th component of $. 

In the next two Propositions we will assume that all the above defined data (the 
Goursat structures T> and T>, the diffeomorphism 0, the coordinates x and x, and 
the Kumpera-Ruiz normal forms K n and £ n ) have been fixed and, therefore, that the 
diffeomorphism $ is uniquely defined. The following result is a direct consequence 
of the obvious relations 

</>*(C») = Cj, 

for < i < n — 4, where C, C X>® denotes the characteristic distribution of 
and Cj that of ( se e Proposition |4.1|) . 



Proposition 7.1 For eac/i 1 < i < 3 we /lave = 3>i(a?i, X2, £3). For each 

4 < i < n we have = $i(a;i, ... , Xj). 

Decompose W 1 into a direct product M n = W x IR™~\ It follows directly from 
Proposition |7.1| that for each 3 < i < n we can build a diffeomorphism 3>W, between 
the projection of x(C/) on R* and the projection of x(?7) on W, by taking the 
components $j, for 1 < j < i, as the components of Denote by the inverse 
of $w. We obviously have $w = $. The following result is a direct consequence 
of Proposition [O and the obvious relations 



which hold for < i < n — 2. Recall that, by definition, the two Kumpera-Ruiz 
normal forms /t" and k n are given by two sequences of prolongations. We will 
denote by n and respectively, the Kumpera-Ruiz normal forms 

obtained as intermediate steps of these successive prolongations. 

Proposition 7.2 There exist four smooth functions, denoted by z/3, 773, /i3 ; and A3, 
that depend on the coordinates X\, X2, and X3 only, such that 

$1 3 )( K 3) = o #(3)^3 + (A 3 o ^( 3 ))£3 
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Moreover, for each i > 4, t/iere exzsi t/iree smooth functions, denoted by Vi, r\i, and 
Hi, that depend on the coordinates x±, . . . ,Xi only, such that 



[Vi o 



fa o + o 



T/ie functions V{, rji, fii, and Aj are uniquely defined, for each i > 3, once the 
diffeomorphism $ has been fixed. They obviously satisfy (v 3 fi 3 — A3ri 3 )(0) ^ and 
(^(0)^0, for i > A. 



rem 



The following result can be considered as a singular version of Backhand's theo- 
see [53] for a modern approach). It shows that any contact transformation is 
the "prolongation" of a first order contact transformation. Though the case n = 4 is 
classical Rl , it seems that our result for n > 5 is new. Notice that a weaker version 



of Theorem [7.3| has already been announced in [ I3J . Independently, an infinitesimal 
version of Theorem |7.3| has been announced in [|7j and proved in |5j] . 



Theorem 7.3 Let be a local (generalized) contact transformation between two 
Goursat structures T> and T>, defined locally at p and p, respectively. Let x and x be 
local coordinates that transform T> and T> into their Kumpera-Ruiz normal forms K n 
and k n , respectively, and let Sx>(p) — w • • ■ u> n -4 be the singularity type of V at p, 
which equals 5f,(p) since V at p and T> at p are locally equivalent. The constants Ci 
and Ci that appear in K n and k n , respectively, and the contact transformation $ 
associated to and to the coordinates x and x fulfill the following relations: 

(i) The diffeomorphism $^ 3 - ) is a first order contact transformation and the func- 
tions v 3 , r] 3 , /i 3 , and A 3 are uniquely determined by 

(ii) The diffeomorphism is uniquely defined by 



$4 {x) 



V3 + ^4??3 
1^3 + Z4A3 



^4 =1^3 + X4X3 

U 4 = L K 4$ 4 = (/i 3 77 3 - A3Z/3)/ (z/ 3 + X 4 7] 3 ) 2 
7/4 = L K 4$ 4 . 
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(iii) Ifi > 5 and Wi-4 7^ a\ then <&W ^ s uniquely defined by 

„ ^_!(0) | Tfc-i(O) 

A«i-i(0) /ii-i(0) 

= ((Xj + Ci)^_i +77i_i) - Ci 

(iv) Ifi>5 and = ai then 3>® ^ s uniquely defined by 



Vi = L re i$i = (/^i-xi/j-i)/ (r/i_! + Xi77i_i) 2 

Therefore, the (generalized) contact transformation $ is uniquely determined by 
the first order contact transformation $( 3 ). 

This Theorem says that any (generalized) contact transformation between two 
Goursat structures is uniquely defined by a first order contact transformation $^ 3 \ 
In fact, the component $4 of is a linear fractional transformation (Mobius 
transformation) whose coefficients are uniquely determined by the components of 
$( 3 ) (compare H]). For i > 5, successively, the component of $W is either, in the 
case of a singular prolongation, a zero preserving linear fractional transformation 
with Xi = being fixed by the fact that 3>W preserves the hypersurface {xi = 0} or, 
in the case of a regular prolongation, by an affine transformation. In both cases, the 
coefficients of the linear fractional transformation or of the affine transformation are 
uniquely determined by 



Proof of Theorem |7.3| If n = 3 then there is nothing to prove. If n = 4 then the 
result is well known (see e.g. 0). Therefore, we can proceed by induction on the 
integer n > 5. Assume that the Theorem is true for n — 1. By Proposition 7/2 , we 
have 



(n-l)/ n -i 
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In other words, the restriction of $ to x(U) H M n_1 , equipped with coordinates 
X\, . . . , x n -\, is a contact transformation between K n_1 and k n ~ l . Since the Theorem 
is assumed to be true for n — 1, each component $j, for 1 < i < n — 1, satisfies the 
relations given by the Theorem, as do, for 3 < % < n — 1, the smooth functions z/j, 



/ii, Aj, and r/i, given by Proposition [f[2| What remains to check is that $ 
and r/ n satisfy our conditions. 

Recall that for any diffeomorphism $ (n) = ($ (n_1) , $ n ) T of R n , such that $( n_1 ) 
depends on the first n — 1 coordinates Xj, . . . , £ n _i only, and for any vector field 
/ = ctf n ~ l + f n on M n , where a is a smooth function on M. n , the vector field / n_1 is 
the lift of a vector field on IR n_1 , and the only non-zero component of f n is the last 
one, we have: 

$!")(/) = (a o *(»))$<»-i)(/»-i) + (( L/ $„) o (45) 

Observe that the vector field ^i™ -1 ^/™" 1 ) is lifted (see Notation |2.1| ) along the 
coordinate x n , which is given by <3> n . 

Regular case: If u>„_4 7^ ai then we have = (x n + Cn)Ki~ + /«2~ ■ This relation, 
together with (|4"5| ) and the induction hypothesis leads to: 

$l")(4) = ((x n + Cn ) o c^-^r 1 ) + 

+ ((L K n$„) O «» 

= (((X„ + C„)^_i + 7M-l) O + (^_! O 

+ ((L K »t)o$ (n) ) K? 

= (/i n _! o *(«)) (( ( "" +c "^;- i 1+,? "- 1 o k™- 1 + 

+ ((L K ,$ n ) o 



By Proposition |7.2| , we know that there exist two smooth functions fi n and r] n (with 



/i n 7^ 0) such that 

$l n) (4) = fan O «? + (/i n o k» 

Comparing the last two relations and taking into account that «™ = g|- while k™ -1 , 
K2 _1 , and ^2 have zeros as components multiplying we see that = L K «$ n . 
From the inductive definition of Kumpera-Ruiz normal forms (regular prolongation) 
given in Section we have 

^2 = {%n "I" C n )ftj "I" ^2 
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We can now conclude that /x n = /i n -i and that 

1 



((x n + c n )i/ n _i + r] 



where 



~ _ ^n-l(O) ?7n-l(0) 



yU n _l(0) /in-l(O) 



Now consider Relation ( fig) gives $l n) ( K «) = ((L„»$ n ) o *(«)) «» which implies 
i/„ = L K n$„. This obviously gives u n = v^xl 



Singular case: If w n _4 = a\ then we have «2 = ^1 1 + L - Together with 



.n-X 



relation (45) and with the induction hypothesis, this relation leads to: 



+ ((L„»$ n ) o *(">) «» 
+ ((L„»$ n ) o «? 



By Proposition |7.2j , we know that there exist two functions // n and ^ n such that 

The same argument as in the regular case implies r\ n = L K n$ n , fx n = v n -\ + x n r] n _i, 
and 



$Jx) 



Xnj-^n—l 



Moreover, like in the regular case, the relation 
implies = L K n$ n = (/^-i^-i)/^! +^j_i) 2 . 



□ 
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7.2 Are Goursat Structures Locally Determined by Their 
Abnormals? 



In this Subsection we will be interested, in the case of Goursat structures, in the 
following question asked by Jakubczyk: "Are nonholonomic distributions deter- 
mined by their abnormal curves?". Several results have been obtained giving a 
positive answer to this question: for stable degenerations of Engel structures by 
Zhitomirskh J73J, for singular contact structures by Jakubczyk and Zhitomirskii |28f . 
for generic distributions of corank at least equal to three, at typical points, by Mont- 
gomery [f|5|]. Recently, Jakubczyk |26| has proved that the answer is positive if we 



consider abnormal curves of the complexified problem, for all distributions with the 
exception of a small subclass. We will show in this Subsection that this subclass 
contains Goursat structures. 

To start with, let us be more precise on what we mean by the statement that 
distributions are determined by their abnormal curves. We will follow the definitions 
given in ||45| . Distributions that belong to a class Q of distributions are strongly 
determined by their abnormal curves if, for any pair of distributions T> and T> that 
belong to Q, any local diffeomorphism that transforms each abnormal curve of T> 
into an abnormal curve of V, and the other way around, transforms also T> into T>. 
It is clear that Goursat structures are not strongly determined by their abnormal 
curves because they have very few abnormal curves. For example, contact structures 
do not have any non-trivial abnormal curve. 

A weaker property can be defined as follows. Distributions that belong to a 
class Q of distributions are weakly determined by their abnormal curves if, for any 
pair of distributions T> and T> that belong to Q, the existence of a local diffeomor- 
phism that transforms each abnormal curve of T> into an abnormal curve of T>, and 
the other way around, implies the local equivalence of T> and T>. 

Proposition 7.4 Goursat structures on n-manifolds are not weakly determined by 
their abnormal curves ifn>6. 



Proof of Proposition |7.4| Consider the two following Kumpera-Ruiz normal forms 
defined on R 6 by 

and 
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On the one hand, by Theorem 1672, the distributions spanned by these two 



Kumpera-Ruiz normal forms have the same abnormal curves, locally at zero. In- 
deed, for each of them, the submanifolds Sj, for j = and 1, are empty in a small 
enough neighborhood of zero (see Proposition (4.4| ); and thus their abnormal curves 
are given, in both cases, by *4.(°) = Cq = (gf-), in a small enough neighborhood of 
zero. But on the other hand, it has been shown by Kumpera and Ruiz |32j that 



these two distributions are not locally equivalent at zero. Indeed, the first one has 
singularity type at zero while the second one has singularity type at 

zero. Analogous examples can be constructed for any n > 6. □ 

Our study of relations between abnormal curves and their singularity type shows 
that the geometry of a Goursat structure is reflected by abnormal curves of all 
elements of the derived flag. It is thus natural to introduce the following definition. 
Distributions that belong to a class Q of distributions are weakly determined by 
abnormal curves of their derived flags if, for any pair of distributions T> and T> that 
belong to Q, the existence of a local diffeomorphism that transforms each abnormal 
curve of T>^ into an abnormal curve of P™, and the other way around, for each 
i > 0, implies the local equivalence of T> and T>. It is a direct consequence of 
Theorem |6.2| and of the classification obtained in |I2]| , ||T9f , and that Goursat 



structures on MP, for 3 < n < 8 are determined by abnormal curves of their derived 
flags. It is surprising that in higher dimensions it is not the case. Indeed, we have 
the following result which is a direct consequence of Theorem |6.2| and the Theorem 



announced in JT3[ 



Proposition 7.5 Goursat structures on n-manifolds are not determined by abnor- 
mal curves of their derived flags if n > 9. 



It has already been announced in p7j| that the growth vector is not a complete 



invariant for Goursat structures on M. n , for n > 9 (which, together with Theorem | 
implies the above result). We will give in this Subsection our proof of this latter 
fact. An alternative proof can be found in J49j. It is important to stress that the 
method used in |49| and the method that we will present in this Section are different. 



It seems that both methods apply, in general, to different cases of non-equivalence. 



Our aim now is to prove Proposition 775. This will be done by giving an example 
(Proposition 7J3) of two Goursat structures T> and T> that are locally non-equivalent 
but that have the same singularity type (which, by Theorem |6.2| , implies the exis- 
tence of a diffeomorphism between the abnormal curves of T>^ and those of P^, 
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for i > 0). Then, this example will be improved (Proposition |7.7| ) by construct- 
ing, instead of a pair of distributions, a continuous family (parametrized by a real 
number) of locally non-equivalent Goursat structures that have the same singularity 
type, and thus diffeomorphic collections of abnormal curves for all elements of their 
derived flags. 

Consider two Kumpera-Ruiz normal forms ( ) and ( ), defined on IR n , 
centered at zero, and given, respectively, in coordinates x = (xi, . . . ,x n ) and 
x = (xi, . . . ,x n ). Assume that they have been obtained from (k\, kQ and (k\, k 1 2 ), 
respectively, by a sequence of regular prolongations, for i > 3. Suppose, moreover, 
that the Goursat structures spanned by (k™, kV,) and (k™, k 2 ) are locally equivalent 
and let x = $(x) be a (generalized) contact transformation, of order n — 2, that 
establishes this equivalence. We have Xj = &j(x), for 1 < j < n. We are going to 
prove that the components $j, for i + 1 < j < n, can be obtained by a sequence of 
derivations (with respect to a well chosen vector field) from the component $j. To 



start with, apply Theorem |7.3| to the component We have 



= — ((Xi+l + Ci + i)Vi + T]i) - Q+i. 
fa 

In follows also from Theorem (regular case) that this expression can be written 
in the following form: 

= — Uxi+i + Q+i)L Kl $i + L K ^i) - Ci+i 
fa v / 

Mi K2 

But since is a function of xi, . . . ,Xi only, the latter expression can be rewritten 

as 

m 2 

Theorem |7.3| implies, moreover, that fa, = fa, for i + 1 < j < n. Thus, the previous 
argument can be repeated to obtain, for 1 < k < n — i, the following relations: 

$ i+fc = lA n $i - c i+k , 

which imply that 



~ Ci+k = iA (o) 
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because the coordinates are centered. Therefore, in the case of a sequence of regu- 
lar prolongations, the constants 5i + k can be obtained by computing the successive 
derivatives L^^.-^^ of the component $j (that defines the coordinate Xj) and by 
taking their values at zero. 

The following definition is natural and will simplify the proofs of the next results 
given in this Subsection. Let 7 be a smooth function defined on R n and let g be a 
smooth vector field, also defined on ~R n . The degree of the function 7, with respect to 
the vector field g, is the smallest integer k (maybe infinite) such that L^(7)(0) 7^ 0. 
Note that if the degree of 71 is i\ and the degree of 72 is ii then the degree of 7172 
is obviously i\ + i 2 - 

Proposition 7.6 Consider the two following Kumpera-Ruiz normal forms defined 
on M 9 by 

^■9 _ _d_ 

1 9xg 

where eg = or 1. They are locally non- equivalent at zero, although both of them 
have the same singularity type CLq CLq CL 1 (I2 ^0 ^0 

at zero. 

Proof of Proposition |7.6| Denote by k 9 the Kumpera-Ruiz normal form given by 
(k\ , ^(0)), in (xi, . . . , ^-coordinates, and denote by k 9 the Kumpera-Ruiz normal 
form given by (k 9 , /^(cg)), in (xi, . . . , Xg)-coordinates. We are going to show that 
if a (generalized) contact transformation x = $(x) converts the Goursat structure 
generated by k 9 into the one generated by k 9 then we must have eg = 0. 

Denote by and by the elements of the two sequences of 

Kumpera-Ruiz normal forms used to construct, via prolongations, the normal forms 
k, 9 and k 9 , respectively. Since k 5 = Rq(k, 4 ) we have, by the regular case of Theo- 
rem |7.tj| , the following relations: 

/i 5 = yU 4 and u 5 = — 
Hence ^5 and z/ 5 are functions of X\,... , X4 only. Denote fi = [14, v = z/ 4 , and 



T] = 775. Since k = S(k ) we have, by the singular case of Theorem |7.3|, the 
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following relations: 

= y 

- + a: 6 ?7 

fi 

v 

He = - + x 6 i]. 
H 

Denote a = l//i 6 and g = olk\. Since both k 9 and k 9 are obtained by a sequence 
of regular prolongations from /t 6 and k 6 , respectively, it follows from the discussion 
given at the beginning of this Subsection that the new constant c g can be calculated 
by computing the successive derivatives of <E> 6 , in the direction of the vector field 
g = (1/// 6 )k|. Namely 

c 9 = (LjS 6 ) (0). 

Instead of computing the successive derivatives of <3> 6 directly, take the Taylor 
series expansion of $ 6 . The terms of this expansion that contain coordinate functions 
of degree d > 4, with respect to g, can obviously be discarded. To this aim, we will 
start by computing the degree, with respect to g, of the functions x±, . . . ,x 6 , that 
is of the variables on which $ 6 depends. 

For xe, we have: 

LgXe = ax-? 

L 2 g x 6 = a 2 (x s + 1) + (L 9 a) x 7 

L 3 g x 6 = a 3 x 9 + 3a (L g a) (x 8 + 1) + (L^a) x 7 . 

Since (L 2 g x 6 ) (0) = a 2 (0) = ( / u(0)/z/(0)) 2 ^ 0, the degree of x 6 is 2. We have 
L g x 5 = a. Therefore the degree of x 5 is 1. We have L s x 4 = ax 6 x 5 . Thus the degree 
of X4 is 4. Analogously, the degree of x 3 is 7, the degree of x 2 is 10, and the degree 
of xi is 3. 

Now observe that $e( x ) = x 6 (p(xi, . . . ,x 6 ), for a suitable function ip. This 
implies that each term of the Taylor series expansion of $ 6 is of the form xex^ 1 ■ ■ ■ x\ 6 , 
for some integers fci, . . . , k e . Since c 9 = (L^$ 6 ) (0), we consider only terms of degree 
d < 3 with respect to g. Therefore we have: 

$e(x) = Axe + Bx 6 x 5 , 

up to terms of degree d > 4 with respect to g. Recall that neither /i nor v depend 
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on the variables X5 and x§. Hence 

<9$ 6 v 



dx 6 

d 2 $ 6 -2x 6 rj xs 



v 

-+X 6 7] 



2 



dx^dxQ 



v 

- + x&ri 



3 • 



Thus A = /i(0) and 5 = 0. This implies that $e(x) = /i(0)x6, up to terms of degree 
d > 4. Since we have already computed the successive derivatives of xq, it is easy 
to obtain that: 

(L,$ 6 ) (0) = 

(L^) (0) = M0)a 2 (0) 

(Lj$ 6 ) (0) = 3/i(0)a(0) (L g a) (0). 

But (L g a) (0) = 0. Hence, since eg = (h^ e ) (0), we have eg = 0. □ 

Proposition 7.7 Consider the following family of Kumpera- Ruiz normal forms de- 
fined on R 11 by 

k 11 = -2- 

1 fen 

where c\\ is an arbitrary real constant. Two Kumpera-Ruiz normal forms that be- 
long to this family are locally equivalent at zero if and only if they have the same 
constant c\\, although all of them have the same singularity type 
at zero. 

Proof of Proposition |7.7| Denote by k 11 the Kumpera-Ruiz normal form given by 
(ki~, ^(cn)), in (xx, . . . , xn)-coordinates, and denote by k 11 the Kumpera-Ruiz 
normal form given by (ftp, /^(cn)), in (£1, . . . , Xn)-coordinates. We are going to 
show that if a (generalized) contact transformation x = $(a;) converts the Goursat 
structure generated by n 11 into the one generated by k 11 then we must have en = en. 

Denote by and by the elements of the two sequences 

of Kumpera-Ruiz normal forms used to construct, via prolongations, the normal 



71 



forms k 11 and re 11 , respectively. Since k 5 = Ro(k, a ) we have, by the regular case of 

•elations: 

//5 = /i4 and z/ 5 



Theorem |7.3| , the following relations: 

Hence /x 5 and z/ 5 are functions of £i,... ,x 4 only. Denote \i = /i 4 , z/ = z/ 4 , and 
= 775. Since = S^k 5 ) we have, by the singular case of Theorem fOj the 
following relations: 

*e(*) = 

V 

/i e = - + x 6 ?7. 
/i 

Denote a = 1///6 an d 9 — OiK^ . Since both k 11 and k 11 are obtained by a sequence 
of regular prolongations from K e and /c 6 , respectively, it follows from the discussion 
given at the beginning of this Subsection that the new constant Cn can be obtained 
by computing the successive derivatives of $6 ; m the direction of the vector field g. 
Namely 

eb = (Lj$ 6 ) (0). 

Let us consider the Taylor series expansion of $6- Again, the terms of this 
expansion that contain coordinate functions of degree d > 6, with respect to g, will 
be discarded. The successive derivatives of Xq are given by: 

IjqXq = axj 

L 2 xq = a 2 x 8 + {L g a)x-j 

L^x 6 = a 3 (x 9 + 1) + 3a(L g a)x 8 + (L 2 a)xj 

L 4 x 6 = a 4 (x 10 + 1) + Qa 2 (L g a)(x 9 + 1) 



+ (3(L 9 a) 2 + 4a(L 2 g a)) x 8 + (L 3 g a)x 7 
L 5 g x 6 = a 5 (x n + en) + 10a 3 (L 9 a)(xi + 1) 
+ (l5a(L 9 a) 2 + 10a 2 (L 2 a)) (x 9 + 1) 
+ (lO(L g a)(L 2 g a) + 5a(L 3 a)) x 8 + (L A g a)x 7 . 

Since (L 3 a; 6 ) (0) = a 3 (0) = (/i(0)/V(0)) 3 7^ 0, the degree of x 6 is 3. The degree of x 5 
is 1, the degree of x^ is 5, the degree of £3 is 9, the degree of £2 is 13, and the degree 
of xi is 4 (all degrees are with respect to g). 
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Now observe that $g(x) = XQip(xi, . . . , xq), for a suitable function ip. The Taylor 
series expansion of $6, up to terms of degree d < 5 with respect to g, is given by: 

$ 6 (x) = Axq + Bxqx 5 + Cxqx\. 

Recall that neither fi nor v depend on the variable X5. Therefore 

<9$ 6 



dx 



: (0) = fi(0) 



6 



^ 6 -(0) = ° (o) =0 . 



(9X5(9X6 <9 2 X5(9X6 

Hence A = /x(0) and both I? and C are equal to 0. This implies that $60*0 = 
/i(0)x 6 , up to terms of degree d > 6. Since we have already computed the successive 
derivatives of xq, it is easy to obtain that: 

(L 9 $ 6 ) (0) = 
(L 2 $ 6 ) (0) = 
(L 3 $ 6 ) (0) = M0)a 3 (0) 
(Lj$ 6 ) (0) = M0)a 4 (0) 
(L*$ 6 ) (0) = //(0)a 5 (0)c n . 

Since $ transforms /t 11 into k 11 , we must have both Cg = 1 and cio = 1. But c 9 = 
(L 3 $ 6 ) (0) and cio = (L^$ 6 ) (0). Therefore, /i 4 (0)/z/ 3 (0) = 1 and / u 5 (0)/z/ 4 (0) = 1. 
This obviously implies /x(0) = z/(0) = 1. Hence, since Cn = Cn/i 6 (0)/z/ 5 (0), we have 
Cn = cn. □ 
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A Weber's Problem 



Our proof of Kumpera-Ruiz's Theorem was based on the following fact: If a rank 
two distribution T> on a manifold M of dimension n > 4 satisfies dimD^Qo) = 3 
and &m\T>( 2 \p) = 4, for each point p in M, then there exists a canonical line field 
£ C T> that satisfies [£, X^ 1 )] C This observation has a natural generalization: 
If a rank k > 2 distribution Dona manifold M of dimension n > k + 2 satisfies 
dim T>^\p) — k + 1 and dimV^ (p) = k + 2, for each point p in M, then there exists 
(i) a canonical involutive distribution C\ C P*- -' that has rank k — 1 and is uniquely 
characterized by [£i,X>W] C P 1 - 1 ^ and (ii) a canonical involutive distribution £ C 
that has rank fc — 2 and is uniquely characterized by [£o ; £^°' ) ] C T>^ (see 152] 



and j£J for an approach based on Pfaffian systems; see also [|3T| and Proposition f4.1| ). 



Though the above observation appears more or less clearly in the work of E. 
Cartan (see e.g. [IT|; see also [23] ), its origin can be found in the pioneering work 



of F. Engel [16], for n = k + 2, and E. von Weber [70], for n > k + 2 (see also |J). 



This observation is closely related to the following result, which is clearly stated 
in Weber's article |70|, Theorem V] (using the dual language of Pfaffian systems). 

Theorem A.l (E. von Weber) Let V be a rank k > 2 distribution on a manifold 
M of dimension n = m + k — 2 > 4. Assume that dim (p) —k + 1 and 
dimX>( 2 )(p) = k + 2, for each point p in M. Then, in a small enough neighborhood 
of any point p in M , the distribution D is equivalent to a distribution spanned by a 
family of vector fields that has the following form: 

m-2 



d d d 



dx m+k _ 2 ' dx m+1 ' dxm ' Xm dx.„^ 

%=1 



where the functions tpi, for 2 < i < m — 2, depend on the variables x±, . . . ,x m -i 
only. 



The following result is a direct consequence of Theorem A.l 



Proposition A. 2 Any Goursat structure on a manifold M of dimension n > 4 is 
equivalent, in a small enough neighborhood of any point p in M, to a distribution 
spanned by a pair of vector fields that has the following form: 

dx~ > Xn dx n . 1 JrXn ~ l dx n ^2 ^ ^^(^n-l)^- + q^I ), (47) 

i=2 / 

where the coordinates centered at p and the functions tpi, for 2 < i < 

n — 3, depend on the variables xi, . . . , x n _i only. 
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In the particular case of four-manifolds the last result gives: 



Corollary A. 3 (Engel's Theorem) Any Goursat structure on a four-manifold 
M is equivalent, in a small enough neighborhood of any point p in M, to the distri- 
bution spanned by the following pair of vector fields (Engel's normal form): 



8 T 8 _i_ „ 8 _i_ 8 
8x4 ' 4 g x:i o q X2 g xi 



where the x-coordinates are centered at p. 

The following Theorem can be considered as a rigorous version of Weber's re- 



sult [ 70 , Theorem VI] . Although it is a direct consequence of the work of Kumpera 



and Ruiz [32], Martin and Rouchon |43|| , and Murray fl5"I |, it seems that it has never 



been stated in the following explicit form. 

Theorem A. 4 (Weber's Problem) A rank k > 2 distribution V on a manifold 
M of dimension n = m + k — 2>4is equivalent, in a small enough neighborhood of 
a given point p in M, to the distribution spanned by the following family of vector 
fields (Weber's normal form) 

8 8 8 „ 8 , „ 8 8 



8x m+k _ 2 ' ' ' ' ' 8x m+1 ' 8x m ' m 8x m - 1 6 8x 2 dx\ 

if and only if dim Vi(p) = dimT>^'(p) = k + i, for < i < m — 2, in a small enough 
neighborhood of p. 

If we have dimP^(p) —k + i, for < i < n — 2, but we do not impose any 
condition on dimDj(p) then we still have the following result, which is a direct 
consequence of Theorem |A.1| and Theorem |2.3| , applied to the last two vectors fields 
of ©. 

Theorem A. 5 (Kumpera-Ruiz) Let V be a rank k > 2 distribution on a mani- 
fold M of dimension n = m + k — 2 > A, such that for any point p in M we have 
dimD^Qo) = k + i, for < i < m — 2. Then, the distribution V is equivalent, in a 
small enough neighborhood of any point p in M , to the distribution spanned by the 
following family of vector fields: 



8x Tl 

where the pair of vector fields (k™, k™) denotes a Kumpera-Ruiz normal form on '. 



75 



B Additional Normal Forms 



Let £ m = be a pair of vector fields defined on M m that has the following 

form: 

m— 3 

^ ' ' t=2 



x ™&z m _i + Xm - 1 ax m _ 2 ~^ S ^(^m-i)^^. + a Xl ) • (49) 



A pair of vector fields £ m+/ = defined on R m+l , for Z > 0, is called a 

prolongation of order I of £ m if we have £ m+l — a\ o ■ ■ ■ o cri(£ m ), where each Cj, for 
1 < i < I, equals either S or R Ci , for some real constants c, (recall that the singular 
and regular prolongations S and R Ci have been defined in Section 0). 
The following Lemma is a natural generalization of Proposition |A.2 . 



Lemma B.l Let V be a rank k > 2 distribution on a manifold M of dimension 
n = m + l + k — 2>4, where I and m are two non-negative integers. Assume that 
for each point p in M we have dimZ>W(p) — ]^ -\- i, for < i < I + 2. Then, in a 
small enough neighborhood of any point p in M , the distribution T> is equivalent to 
a distribution spanned by a family of vector fields that has the following form: 



tin+i+l 

where the pair of vector fields , is a prolongation of order I of a pair of 

vector fields (£™, £™) of the form 



The proof of Lemma |B.1| is left to the reader. For generic points, the Lemma is 



stated and proved in the work of Cartan [|ll[| and Goursat ||22|| . For singular points, 



the Lemma is a direct consequence of the results obtained by Kumpera and Ruiz |32[ 
and its proof is almost the same as that of Theorem |2]3 but there are essentially 



two differences. The first difference is that instead of using Proposition 2.5, as it is 



done in the Proof of Theorem one uses Theorem [A.l| ; the second difference is 



that instead of starting the induction argument, for / = 0, with the Pfaff-Darboux 



normal form, as it is done in the Proof of Theorem [2.3| , one starts it with Weber's 



preliminary normal form 

Let T> be a rank k > 2 distribution on a manifold M of dimension n = m+k—2 > 
4, such that for any point p in M we have diml)W(p) = k + i, for < i < m — 2. 
It is easy to check that each distribution T>^\ for < i < m — 4, contains a 
unique involutive subdistribution Q C T>^ that has constant corank one in 
and is characteristic for D^ +1 \ We can generalize the canonical submanifolds Sq^ 
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of Section |] by the following definition: 

Sit ) = {peM:V^(p)=C l+l (p)}, 

where < % < m — 5. We say that a point p of M is singular if there exists some 
< i < m — 5 such that p G Sq\ For a singular point p, we denote by ko the 
smallest integer 1 < k < m — 4 such that p G S^ ^. 

Lemma B.2 Let V be a rank k > 2 distribution on a manifold M of dimension 
n = m + ko + k — 2>A, such that for any point p in M we have dim (p) = k + i, 
for < i < k + 2. Assume, moreover, that ko is the smallest integer such that 

in a small enough neighborhood ofp, the distribution D is 
equivalent to a distribution spanned by a family of vector fields that has the following 
form: 

( 9 8 pm+ko £vn+k \ 

\ dx m+ko+k _ 2 ' ' ' ' ' dx m+ko+1 'SI 5 <52 J ' 

where the pair of vector fields £ m+k ° = a prolongation of order 

ko of a pair of vector fields £ m = (^] Ti ,^™) of the form (|49|). Moreover, we have 
^m+ko _ a ^ o ■ • ■ o (Ti(^ m ), where a 1 = S and each aj, for 2 < j < k , equals R Cj , 
for some real constants Ci. 

The proof of Lemma |B . 2| follows the same line as that of Proposition [4.4| . Though 
instead of considering a Kumpera-Ruiz normal form we consider now a family of 
vector fields of the form ([51]), the idea is the same. Firstly, we compute the distribu- 
tions a nd Ci, and the submanifolds Sq\ Secondly, we observe that if o\ = R c , 
for some real constant c, then p ^ S^ 0-1 - 1 ; since p G S^ ' 1 ^ we must have <7i = S. 
Thirdly, we observe that if aj = S for some 2 < j < k then p G S^ , since /c is 
by definition the smallest integer such that p G So we must have = R Cj , for 
2 < J < fco- 

For Goursat structures, using the singularity type leads to the following stronger 
result, which states that if the singularity type is of the form waic^ ■ • ■ ctfc then the 
constants that appear in all regular prolongations in the above Lemma equal zero. 

Lemma B.3 Let V be a Goursat structure on a manifold M of dimension n > 5 
andletp beapointinM. If the singularity type of T> atp is of the formwa\a2 - ■ ■ a^, 
for some 1 < k < n — 4, where w is an arbitrary word of J n -fc -3? then T> is locally 
equivalent to a distribution spanned by a pair of vector fields that has the following 
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form: 



9xi 



V i=fco+5 / 

where the coordinates Xi, . . . ,x n are centered at p. 

The proof of the last Lemma follows also the same line as the Proof of Proposi- 
tion |4.4j. Again, we leave details to the reader. The main interest of the last Lemma 



is that it gives directly the proof of Lemma 5.9 



Proof of Lemma |6.9| It is straightforward to check that, in the coordinates of 

-,(feo-l 
'fco-l 



Lemma [B.3|, the canonical sub manifold Si k °_^ is given by 



S^ = {x 1 = 1 ...,x ho = 0} 

and that, moreover, we have Co = (£i) on M and A^^ip) = (^)(p) for each point p 
on Sfo~i- In order to obtain the required normal form, we only have to change 
two coordinates. For 1 < i < k + 2 and k + 5 < % < n, take yi = Xj. Moreover, 
take |/fc +4 — an d yk +3 — — ^fco+s^fco+i + 2 X ka+2 x k +i- n 
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C Figures of Low- Dimensional Trailer Systems 
C.l The Unicycle and the Car 




X\ = U\ 
x 2 = xi u 2 
X3 = U 2 



Figure 1: The unicycle and its normal form. Growth vector: (2,3). Singularity 
type: e. 




Xi 

x 4 



Ui 

Xi u 2 
x 2 u 2 
u 2 



Figure 2: The car and its normal form. Growth vector: (2,3,4). Singularity type: a . 
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C.2 The Two- Trailer System 




X\ 

x 2 
x 3 
X4 
x 5 



Ui 

x 1 u 2 
x 2 u 2 
X3 u 2 
u 2 



Figure 3: A two-trailer and its normal form. Growth vector: (2,3,4,5). Singularity 
type: a a . 




Xi 

x 2 

X3 
X4 
x 5 



Ui 

u 2 

Xi x 2 u 2 
Xi x 3 u 2 
x 1 u 2 



Figure 4: A two-trailer and its normal form. Growth vector: (2,3,4,4,5). Singu- 
larity type: a ai. 
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C.3 The Three- Trailer System 



X\ = U\ 

x 2 = X 1 u 2 

x 3 = x 2 u 2 

x 4 = x 3 u 2 

x 5 = X A u 2 

X 6 = U 2 

Figure 5: A two-trailer and its normal form. Growth vector: (2,3,4,5,6). Singu- 
larity type: a a a . 





Xi 

x 2 
X4 

X 5 

x 6 



Ui 

u 2 

x 1 x 2 u 2 
Xi x 3 u 2 
X\ X A u 2 
X 1 u 2 



Figure 6: A two-trailer and its normal form. Growth vector: (2,3,4,4,5,5,6). 
Singularity type: a^aoai. 
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Figure 7: A two-trailer and its normal 
larity type: a aia . 



X\ = U\ 

£2 = (x 1 + l)u 2 

x 3 = u 2 

X 4 = X 2 X3 U 2 

£5 = x 2 x 4 u 2 

x 6 = X 2 u 2 

. Growth vector: (2,3,4,5,5,6). Singu- 




Xi 

x 2 
x 3 

Xi 

x 5 

Xq 



Xi u 2 
u 2 

x 2 x 3 u 2 
x 2 X A u 2 
x 2 u 2 



Figure 8: A two-trailer and its normal form. Growth vector: (2,3,4,5,5,5,6). 
Singularity type: a aia 2 . 
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X\ = U\ 

X 2 = U 2 

X 3 = Xi U 2 

X4 = X\ X2 X3 U2 

X 5 = Xi X 2 X4 U 2 

Xq = X 1 X 2 U 2 

Figure 9: A two-trailer and its normal form. Growth vector: (2,3,4,4,5,5,5,6). 
Singularity type: 




83 



References 



[1] A. Agrachev and A. Sarychev. On abnormal extremals for Lagrange variational 
problems. Journal of Mathematical Systems, Estimation, and Control, 8(1):87- 
118, 1998. 

[2] A. Backlund. Uber Flachentransformationen. Mathematische Annalen, 9:297- 
320, 1876. 

[3] B. Bonnard and I. Kupka. Theorie des singularites de l'application entree/sortie 
et optimalite des trajectoires singulieres dans le probleme du temps minimal. 
Forum Mathematicum, (5):111-159, 1993. 

[4] R. Brockett. Control theory and singular Riemannian geometry. In P. Hilton 
and G. Young, editors, New Directions in Applied Mathematics, pages 11-27. 
Springer- Verlag, New York, 1981. 

[5] R. Brockett. Asymptotic stability and feedback stabilization. In R. Brockett, 
R. Millman, and H. Sussmann, editors, Differential Geometric Control Theory, 
pages 181-191. Birkhauser, Boston, 1983. 

[6] R. Bryant, S-S. Chern, R. Gardner, H. Goldschmidt, and P. Griffiths. Exterior 
Differential Systems. Mathematical Sciences Research Institute Publications. 
Springer- Verlag, New York, 1991. 

[7] R. Bryant and L. Hsu. Rigidity of integral curves of rank 2 distributions. 

Inventiones Mathematicae, (114):435-461, 1993. 

[8] M. Cahadas-Pinedo and C. Ruiz. Pfafhan systems with derived length one. 
The class of flag systems. Preprint, University of Granada (Spain). 

[9] E. Cartan. Sur l'integration de certains systemes de Pfaff de caractere deux. 
Bulletin de la Societe Mathematique de France, 29:233-302, 1901. (Euvres 
completes, Part. II, Vol. 1, Gauthiers-Villars, Paris. 

[10] E. Cartan. Les systemes de Pfaff a cinq variables et les equations aux derivees 
partielles du second ordre. Annales de I'Ecole Normale Superieure, 27:108-192, 
1910. Giuvres completes, Part. II, Vol. 2, Gauthiers-Villars, Paris. 

[11] E. Cartan. Sur l'equivalence absolue de certains systemes d'equations 
differentielles et sur certaines families de courbes. Bulletin de la Societe 

84 



Mathematique de France, 42:12-48, 1914. (Euvres completes, Part. II, Vol. 
2, Gauthiers-Villars, Paris. 

[12] M. Cheaito and P. Mormul. Rank-2 distributions satisfying the Goursat con- 
dition: All their local models in dimension 7 and 8. ESAIM Control, Optimi- 
sation, and Calculus of Variations, 4:137-158, 1999. 

[13] M. Cheaito, P. Mormul, W. Pasillas-Lepine, and W. Respondek. On local clas- 
sification of Goursat structures. Comptes Rendus de VAcademie des Sciences 
de Paris, 327(5) :503-508, 1998. 

[14] J.-M. Coron. Global asymptotic stabilization for controllable systems without 
drift. Mathematics of Control, Signals, and Systems, 5:295-312, 1991. 

[15] G. Darboux. Sur le probleme de Pfaff. Bulletin des Sciences mathematiques, 
2(6):14-36,49-68, 1882. 

[16] F. Engel. Zur Invariantentheorie der Systeme Pfaff 'scher Gleichungen. Berichte 
iiber die Verhandlungen der Koniglich Sdchsischen Gesellshaft der Wissen- 
shaften. Mathematisch-Physikalische Klasse, Leipzig, 41,42:157-1 76;192-207, 
1889,1890. 

[17] M. Fliess, J. Levine, P. Martin, and P. Rouchon. Flatness and defect of non- 
linear systems: Introductory theory and examples. International Journal of 
Control, 61(6):1327-1361, 1995. 

[18] G. Frobenius. Uber das Pfaff 'sche problem. Journal fur die reine und ange- 
wandte MathemaUk, 82:230-315, 1877. 

[19] M. Gaspar. Sobre la clasificacion de sistemas de Pfaff en bandera. In Proceedings 
of the Spanish-Portuguese Conference on Mathematics, pages 67-74, Murcia 
(Spain), 1985. 

[20] A. Giaro, A. Kumpera, and C. Ruiz. Sur la lecture correcte d'un resultat d'Elie 
Cartan. Comptes Rendus de VAcademie des Sciences de Paris, 287:241-244, 
1978. 

[21] E. Goursat. Sur le probleme de Monge. Bulletin de la Societe Mathematique 
de France, (33):201-210, 1905. 

[22] E. Goursat. Legons sur le probleme de Pfaff. Hermann, Paris, 1923. 



85 



[23] D. Hilbert. Uber den Begriff der Klasse von Differentialgleichungen. Mathema- 
Usche Annalen, 73:95-108, 1912. 

[24] B. Jacquard. Le probleme de la voiture a deux, trois et quatre remorques. 
Preprint, DMI-ENS Paris, 1993. 

[25] B. Jakubczyk. Invariants of dynamic feedback and free systems. In Proceedings 
of the European Control Conference, pages 1510-1513, Groningen, 1993. 

[26] B. Jakubczyk. Characteristic varieties of distributions and abnormal curves. 
Preprint, 1999. 

[27] B. Jakubczyk and F. Przytycki. Singularities of /c-tuples of vector fields. Dis- 
sertationes Mathematicae, (213): 1-64, 1984. 

[28] B. Jakubczyk and M. Zhitomirskh. Odd-dimensional Pfaffian equations: Re- 
duction to the hypersurface of singular points. Comptes Rendus de I'Academie 
des Sciences de Pans, 325(4) :423-428, 1997. 

[29] F. Jean. The car with n trailers: Characterization of the singular configurations. 
ESAIM Control, Optimisation, and Calculus of Variations, 1:241-266, 1996. 

[30] Z-P. Jiang and H. Nijmeijer. A recursive technique for tracking control of non- 
holonomic systems in chained form. IEEE Transactions on Automatic Control, 
44(2):265-279, 1999. 

[31] M. Kazarian, R. Montgomery, and B. Shapiro. Characteristic classes for the 
degenerations of two-plane fields in four dimensions. Pacific Journal of Math- 
ematics, 179(2):355-370, 1997. 

[32] A. Kumpera and C. Ruiz. Sur l'equivalence locale des systemes de Pfaff en 
drapeau. In F. Gherardelli, editor, Monge- Ampere equations and related topics, 
pages 201-247. Instituto Nazionale di Alta Matematica Francesco Severi, Rome, 
1982. 

[33] G. Lafferriere and H. Sussmann. A differential geometric approach to motion 
planning. In Z. Li and J. F. Canny, editors, Nonholonomic motion planning, 
International Series in Engineering and Computer Sciences, pages 235-270. 
Kluwer, Dordrecht, 1992. 



86 



[34] J. -P. Laumond. Controllability of a multibody mobile robot. IEEE Transactions 
on Robotics and Automation, 9(6):755-763, 1991. 

[35] J. -P. Laumond. Singularities and topological aspects in nonholonomic motion 
planning. In Z. Li and J. F. Canny, editors, Nonholonomic motion planning, 
International Series in Engineering and Computer Sciences, pages 755-763. 
Kluwer, Dordrecht, 1992. 

[36] J. -P. Laumond, editor. Robot Motion Planning and Control. Lecture Notes on 
Control and Information Sciences. Springer- Verlag, Berlin, 1997. 

[37] J.-P. Laumond, P. Jacobs, M. Taix, and R. Murray. A motion planner for 
nonholonomic mobile robots. IEEE Transactions on Robotics and Automation, 
10(5):577-593, 1994. 

[38] Z. Li and J.-F. Canny, editors. Nonholonomic Motion Planning. International 
Series in Engineering and Computer Sciences. Kluwer, Dordrecht, 1992. 

[39] P. Libermann. Sur le probleme d'equivalence des systemes de Pfaff non 
completement integrables. Publications Paris VII, 3:73-110, 1977. 

[40] S. Lie and G. Scheffers. Geometric of Beriihrungstransformationen. B. G. 
Teubners, Leipzig, 1896. 

[41] W. Liu. An approximation algorithm for non-holonomic systems. SIAM Journal 
on Control and Optimization, 35(4):1328-1365, 1997. 

[42] F. Luca and J-J. Risler. The maximum degree of nonholonomy for the car 
with n trailers. In Proceedings of the IFAC Symposium on Robot Control, pages 
165-170, Capri (Italy), 1994. 

[43] P. Martin and P. Rouchon. Feedback linearization and driftless systems. Math- 
ematics of Control, Signals, and Systems, 7:235-254, 1994. 

[44] R. M'Closkey and R. Murray. Exponential stabilization of driftless nonlinear 
control systems using homogeneous feedback. IEEE Transactions on Automatic 
Control, 42(5):614-628, 1997. 

[45] R. Montgomery. A survey of singular curves in sub-Riemannian geometry. 
Journal of Dynamical and Control Systems, (l):49-90, 1995. 



87 



[46] P. Morin and C. Samson. Exponential stabilization of nonlinear driftless sys- 
tems with robustness to unmodeled dynamics. ESAIM Control, Optimisation, 
and Calculus of Variations, 4:1-35, 1999. 

[47] P. Mormul. Contact hamiltonians distinguishing locally certain Goursat sys- 
tems. Preprint, Warsaw, 1998. 

[48] P. Mormul. Local models of 2-distributions in 5 dimensions everywhere fulfilling 
the Goursat condition. Research report, Rouen, 1994. 

[49] P. Mormul. Rank-2 distributions satisfying the Goursat condition: All their lo- 
cal models in dimension 9. Preprint, Institute of Mathematics - Polish Academy 
of Sciences, 1997. 

[50] P. Mormul. Goursat distributions with one singular hypersurface - constants 
important in their Kumpera-Ruiz pseudo- normal forms. Preprint, Universite 
de Bourgogne, 1999. 

[51] R. Murray. Nilpotent bases for a class of nonintegrable distributions with 
applications to trajectory generation for nonholonomic systems. Mathematics 
of Control, Signals, and Systems, 7:58-75, 1994. 

[52] R. Murray and S. Sastry. Nonholonomic motion planning: Steering using sinu- 
soids. IEEE Transactions on Automatic Control, 38(5):700-716, 1993. 

[53] P. Olver. Equivalence, Invariants, and Symmetry. Cambridge University Press, 
1995. 

[54] W. Pasillas-Lepine and W. Respondek. Applications of the geometry of Gour- 
sat structures to nonholonomic control systems. In Proceedings of the IF AC 
Nonlinear Control Systems Design Symposium, pages 789-794, Enschede (The 
Netherlands), 1998. 

[55] W. Pasillas-Lepine and W. Respondek. Conversion of the n-trailer into 
Kumpera-Ruiz normal form and motion planning through the singular locus. 
In Proceedings of the IEEE Conference on Decision and Control, Phoenix (Ari- 
zona), 1999. To appear. 

[56] J.-B. Pomet. Explicit design of time- varying stabilizing control laws for a class 
of controllable systems without drift. Systems and Control Letters, 18:147-158, 
1992. 



88 



[57] L. Pontryagin, V. Boltyanskff, R. Gamkrelidze, and E. Mischenko. The Math- 
ematical Theory of Optimal Processes. Wiley, New York, 1962. 

[58] P. Rouchon, M. Fliess, J. Levine, and P. Martin. Flatness and motion planning: 
The car with n trailers. In Proceedings of the European Control Conference, 
pages 1518-1522, Groningen, 1993. 

[59] C. Samson. Control of chained systems: Application to path following and time- 
varying point-stabilization of mobile robots. IEEE Transactions on Automatic 
Control, 40(l):64-77, 1995. 

[60] O. S0rdalen. Conversion of the kinematics of a car with n trailers into a chained 
form. In Proceedings of the IEEE Conference on Robotics and Automation, 
pages 382-387, Atlanta (Georgia), 1993. 

[61] O. S0rdalen. On the global degree of nonholonomy of a car with n trailers. In 

Proceedings of the IFAC Symposium on Robot Control, pages 343-348, Capri 
(Italy), 1994. 

[62] O. S0rdalen and O. Egeland. Exponential stabilization of nonholonomic chained 
systems. IEEE Transactions on Automatic Control, 40(l):35-49, 1995. 

[63] O. S0rdalen, Y. Nakamura, and W. Chung. Design and control of a non- 
holonomic manipulator. In Ecole d'ete d'automatique de VENSIEG, Grenoble 
(France), 1996. 

[64] O. S0rdalen and K. Wichlund. Exponential stabilization of a car with n trailers. 

In Proceedings of the IEEE Conference on Decision and Control, pages 978-983, 
San Antonio (Texas), 1993. 

[65] H. Sussmann and W. Liu. Shortest paths for sub-Riemannian metrics of rank-2 
distributions. Memoirs of the American Mathematical Society, 192(118), 1995. 

[66] A. Teel, R. Murray, and G. Walsh. Nonholonomic control systems: From steer- 
ing to stabilization with sinusoids. International Journal of Control, 62(4) :849- 
870, 1995. 

[67] D. Tilbury, R. Murray, and S. Sastry. Trajectory generation for the n-trailer 
problem using Goursat normal form. IEEE Transactions on Automatic Control, 
40(5):802-819, 1995. 



89 



[68] A. Vershik and V. Gershkovich. Nonholonomic dynamical systems, geometry of 
distributions and variational problems. In V. Arnol'd and S. Novikov, editors, 
Dynamical systems VII, Encyclopaedia of Mathematical Sciences. Springer- 
Verlag, New- York, 1991. 

[69] G. Walsh, D. Tilbury, S. Sastry, R. Murray, and J-P. Laumond. Stabilization 
of trajectories for systems with nonholonomic constraints. IEEE Transactions 
on Automatic Control, 39(l):216-222, 1994. 

[70] E. von Weber. Zur Invariantentheorie der Systeme Pfaff'scher Gleichungen. 

Berichte iiber die Verhandlungen der Koniglich Sachsischen Gesellshaft der 
Wissenshaften. Mathematisch-Physikalische Klasse, Leipzig, 50:207-229, 1898. 

[71] I. Zelenko and M. Zhitomirskh. Rigid paths of generic 2-distributions on 3- 
manifolds. Duke Mathematics Journal, 79(2):281-307, 1995. 

[72] P. Zervos. Le probleme de Monge. Memorial des Sciences Mathematiques. 
Gauthier-Villars, Paris, 1932. 

[73] M. Zhitomirskh. Normal forms of germs of 2-dimensional distributions on R 4 . 
Functional Analysis and Applications, 24(2): 150-152, 1990. 

[74] M. Zhitomirskh. Normal forms of germs of distributions with a fixed segment of 
growth vector. Leningrad Mathematics Journal, (2):1043-1065, 1991. (English 
translation). 

[75] M. Zhitomirskh. Rigid and abnormal line subdistributions of 2-distributions. 

Journal of Dynamical Control Systems, (l):253-294, 1995. 

[76] M. Zhitomirskh. Singularities and normal forms of smooth distributions. In 
B. Jakubczyk, W. Respondek, and T. Rzezuchowski, editors, Geometry in Non- 
linear Control and Differential Inclusions, pages 395-409. Banach Center Pub- 
lications, Warsaw, 1995. 



90 



